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SUMMARY 


ANALYTICAL EVALUATION OF ILM SENSORS 
(Supplement) 

Raymond J. Kirk 
C, Ed Johnson 
Douglas Doty 
Honeywell Inc . 

This report is a supplement to CR-132687 , Analytical Evaluation of ILM 
Sensors, September 1975. It documents the analysis and evaluation 
of the capabilities of three specific sensing concepts to meet the 
necessary functional requirements and operational environment for 
an Independent Landing Ifonitor (ILM) which will provide the flight 
crew with an independent assessment of the operation of the primary 
landing system during Category II/III operations. The three concepts 
are; an airborne bistatic radar, a complex interferometer/DME combin- 
ation and a rotating synthetic aperture system. Functional require- 
ments and environmental conditions for the ILM application are 
documented in CR-132687. 

Using our baseline functional requirements, the most critical need for 
the ILM appears in the landing sequence at 1000 to 2000 meters from 
threshold (middle marker) through rollout. Each of the concepts ^ 

evaluated in this report show potential of becoming feasible ILM con- 
cepts. The complex interferometer/DME concept has been developed as 
a Cat. II landing aide and with increases to the ground equipment 
conqjlex should also be capable of some Cat. Ill ILM functions. The other 
two concepts are in early development stages at this time. 
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SECTION I 
INTRODUCTION 


The technology required for automatic landing of conventional fixed-wing air- 
craft using the VHF Instrument Landing System (ILS) and highly redundant flight 
control systems in Category III Conditions has been demonstrated many times. How 
ever, the introduction of this capability into general usage has been limited. 

One reason for the limited use of autoland systems has been the lack of a mon- 
itoring system operating independently from the primary guidance system. Inde- 
pendent operation would allow the pilot to assess the performance of the autoland 
system and the aircraft situation relative to the runway under conditions of low 
visibility and/or allow the pilot to manually land the aircraft in the event of 
failures in a fail-passive autoland system. This independent monitoring function 
should serve to build up the pilot's confidence in low visibility landings and 
permit the reduction of the flight control system from a fall-operational auto- 
matic unit to a fall-passive system with manual backup capability. The Indepen- 
dent Landing Monitor (ILM) is here defined to include the sensors, signal pro- 
cessors, and displays needed to provide the pilot with information to assess his 
situation and either make a missed approach or proceed in with thenanual landing. 

The primary objective of the ILM study was to evaluate analytically the capabil- 
ities of various generic sensor concepts to meet the functional requirements for 
an Independent Landing Monitor. 

To evaluate the capabilities of the various sensing concepts over a range of ILM 
requirements, the functions required of the ILM were defined for a basic con- 
figuration and for three configurations of increasing complexity: 

a. Basic - The ILM shall provide information to allow determination of aircraft 
position (elevation, azimuth, range) with respect to the runway down to a 
wheel altitude of 12 feet. The accuracy of the position Information should 
be adequate to monitor the approach with a Category III Microwave Landing 
System (MLS) as the primary guidance system. 
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b. Case I - The requirements for the basic system configuration are extended to 
touchdown and through rollout to provide guidance to the aircraft through 
the full landing sequence. Heading with respect to end of runway should be 
provided with an accuracy of 0.1°. 

c. Case II - Add the detection of obstacles, such as aircraft and ground vehicles 
on or near the runway to the Basic and Case I configuration functions. 

Detection of obstacles should be at a range sufficient to permit execution 
of a missed approach. 

d. Case III - Meet the Basic, Case I, Case II requirements through the use of a 
real-world perspective image of the runway and its surroundings. The image 
should be of sufficient quality to provide position, rate of change of posi- 
tion, and attitude information to the pilot in a form analogous to the visual 
cues available during VFR approach and landing. 

This report documents an add-on analysis effort to the ILM study which was initiated 
after compilation of the final report was completed. The add-on effort consisted 
of analysis and evaluation of three specific sensing concepts for ILM not included 
in the primary study and summation of the results. 

The final report for the primary study is titled "Analytical Evaluation of ILM 
Sensors", CR-132687, Sept. 1975, It contains a description of the study technical 
approach, summary of results, detail description of ILM functional requirements 
and operational environment, analysis and evaluation on twenty potential ILM 
sensing concepts and further refinement on six concepts considered feasible. 

This add-on effort utilized the technical approach, functional requirements and 
operational environment and some of the computer programs and data already docu- 
mented in CR-132687. This is not duplicated here and readers should refer to 
that report for this information. References to specific sections of that report 
are contained where applicable in this document. 
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This report contains two major sections for ease of reading. A summary of results 
section provides a brief description of the three ILM sensing concepts reviewed, 
a summary of their capabilities and constraints and conclusions and recommendations. 
The second section contains the documentation for the detail analysis and results 
for each of the three concepts. 



SECTION II 
SUMMARY OF RESULTS 


Landing safety and reliability, independent of terminal area weather conditions, 
is a goal towards which the air transport industry is striving. The benefits to 
passengers, the airlines and the military are clear. To this end, fail -opera- 
tional landing systems are minimum requirements for achieving all-weather landing 
Capabilities. The systems, utilizing an accurate and reliable microwave landing 
system as the primary landing aid, form the basis for the FAA's and ICAO's plans 
to make Category III operations available at a limited number of high density 
airports in the 1980's. 

To date, fail-operational systems for Category III conditions have required 
highly-redundant fail-operational flight control systems and high reliability 
landing aids in fully automatic configurations. An alternative approach to the 

■ ■ i 

fully automatic configurations is to utilize a fail-passive flight control system 
and primary landing aid with an independent back-up landing aid and manual flight 
control. The fail-passive primary landing system permits the use of dual re- 
dundancy rather than the more costly triple or quad redundancy of an automatic 
fail-operational system. The back-up independent landing monitor (ILM) with 
manual takeover capability upgrades the system from fall- passive to fall-operational. 

The intent of the Analytical Evaluation of ILM Sensors study was to survey potent- 
ial sensors and sensing concepts to determine their utility as the independent 
landing aid and to identify the limits or conditions under which they would ful- 
fill one or more of the required functions. Development of a baseline set of 
functional requirements and environmental operating conditions for an ILM and 
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evaluation of the feasibility of twenty different sensing concepts to meet these 
requirements was undertaken in the study and is documented in Reference 1. 

The intent of this add-on effort was to analyze and evaluate against the same 
baseline criteria three additional specific sensing concepts. Each of these 
concepts appeared on preliminary examination to be potentially feasible tlM 
candidates, but were identified after completion of the technical effort on the 
primary study. Evaluation of the three concepts to meet the criteria of Reference 
1 is reported here. 

The three concepts are an airborne bistatic thinned array radar under development 
by Hughes Aircraft, Inc. with the tradename "BISTAR", a complex interferometer/ 

DME system in late development as a Category II landing aide by M.E.L., Ltd. with 
the tradename "MADGE", and a rotating antenna synthetic circular aperture device 
under development at ONERA with the tradename "Astrolabe". 

AIRBORNE BISTATIC THINNED ARRAY RADAR 
Concept Description 

This airborne ILM concept uses bistatic radar techniques to achieve improvements 
in resolution. Airborne radars commonly. use the same antenna for both transmit 
and receive. In this case, the transmit antenna is mounted in the nose of the 
aircraft and scans through the required coverage with a relatively broad beam an- 
tenna pattern. A thinned array of receiving elements spread across the wing 
generates a multi-lobe receiving pattern with each lobe having high azimuth 
resolution. The superposition of the scanned transmit antenna pattern and the 
receiving pattern results in a high resolution azimuth scan. Because a separate 
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transmit antenna is used, the wing elements can be separated by approximately 
one-half of the transmit antenna diameter rather than a fraction of the wave- 
length. This results in on the order of 20-30 receiving elements on the wing 
rather than hundreds of transmit-receive elements as would be required for an 
equivalent monostatic antenna array. 

A sketch of the concept is shown in Figure 1. The system is an air- derived 
concept and could be a completely airborne system although passive ground augmen- 
tation to enhance the runway outline may be desirable. The frequency of operation 
could be X, Ku, or Ka band. A key feature of this concept is the ability to 
serve other functions, such as weather detection, in addition to the ILM function. 
Thus, frequency selection may depend also on these other modes. 

The simplest display for use with the bistatic thinned array radar concept is a 
sector PPI image using range-azimuth format. Techniques to generate a perspective 
image with this concept should be explored. 

Summary of Capabilities 

The analysis showed this application has potential for providing high resolution 
direct imaging of the runway during approach in severe rain, snow, or zero visi- 
bility fog. Direct imaging of the runway is limited to within 5Km even for high 
contrast runway detection without enhancement. Specific results are itemized below: 

1. Lowest backscattering cross-sections provide inadequate background to runway 

backscattering cross-section ratio contrast, AB, and background to precipitation 
clutter, B /C, to detect the runway unless circular polarization is used for 
clutter reduction and degraded range resolution is used to decrease the noise 
bandwidth. With these assumptions, approximate detection ranges for 20 db 
contrast (runway to background backscattering cross-section ratio) are: 


Weather 

Detection Range, 23m runway 

16mm/ hr Rain 

1 .9 km 

9mm/ hr Snow 

2.0 km 

Fog, zero visibility 

1.9 km 
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AIRBORNE 


GROUND 



CORNER REFLECTORS 
MAY BE USED TO 
ENHANCE RUNWAY 
OUTLINE 


• FREQUENCY X, Ku, Ka 

• SENSOR AND CONTROL AIRBORNE 


• MEASUREMENT IMAGE 


• SIMPLEST DISPLAY SECTOR PPI 


Figure 1. Airborne Bistatic Thinned Array Radar 






2. Runway images derived by passive reflector enhancement or imaging of cultural 
objects such as runway lights or other objects with approximately one square 
meter scattering cross-sections can be detected with the extended detection 
ranges shown: 


Weather 

Detection Range 

16mm/hr Rain 

2.6 km 

9mm/hr Snow 

3.0 km 

zero visibility 

2.8 km 


3. Higher backscattering cross-sections allow further utilization of the high 
resolution capabilities of the Marconi illuminator with the array receiving 
System. For this case, a runway- to- background backscattering cross-section 
ratio contrast of 35 dB was analyzed. The approximate detection ranges are: 


Weather 

Detection Range, 23m runway 

16mm/ hr Rain 

4.1 km 

9mm/hr Snow 

5.4 km 

Fog, zero visibility 

4.4 km 


4. The limitations above for the 23 meter runway are due to the runway to background 

contrast. This contrast of 35 dB is approximately the maximum allowable based 

on the 33 dB sidelobes measured for the array (Reference 2). Using this as 

the limiting case, assuming 6 dB B /C minimum required for detection, the 

S 

niaximum precipitation rates may be derived for the detection of a 23 meter 
runway at 1 nm and 2 nm. These maximum rates are derived from the approximate 
expression Clutter Power, C ® R^'^, where R is the precipitation fall rate 
in mm/hr. (Reference 1) 
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Maximum Precipitation Rates for 23 meter Runway Detection 


INM 


2nM 


Rain 67 mm/hr 

Snow 59 nun/ hr 

Rain 20 mm/hr 

Snow 17 mm/ hr Equivalent Water 


5. For comparison, the Ka band perspective imaging radars examined in Reference 1 
show detection ranges of 1 to 1.5 km and 1 to 2 km for the conditions of item- 
ized results 1 and 2, respectively. Thus, nearly a factor of 2 Improvement in 
range is shown. However, it is still well below desired ranges. 

COMPLEX INTERFEROMETER/DME 
Concept Description 

This concept uses interferometric techniques to generate elevation /azimuth angle 
data and DME techniques for range measurements in the landing sequence. A sketch 
and block diagram of the concept is shown in Figure 2. The aircraft equipment 
consists of an interrogator and a receiver-decoder. Interrogations are initiated 
by the aircraft system to the ground transponder. The interrogation message is 
received by the transponder and the elevation and azimuth angle antenna units 
On the ground. The elevation and azimuth units utilize complex space diversity 
interferometer arrays to accurately measure angles to the interrogating aircraft* 
in the presence of multi-path signals. The angle data is transferred to the 
transponder and transmitted back to the aircraft as part of the response. Range 
is measured on the aircraft by. measuring the time from interrogation to response. 

Frequency of operation used is C-band. Measurements consist of range- to -touchdown 
and azimuth/ elevation angles. The simplest display will be steering commands and 
range-to-go to touchdown. 

The particular application analyzed here is the MADGE system, developed by the MEL 
Equipment Co. Ltd, Sussex, England. It utilizes a space diversity interferometer to 
obtain elevation data, a standard interferometer with a special array ambiguity 
resolution technique to obtain azimuth, and a pulse DME to obtain range. The pre- 
vious analysis of Independent Landing Monitor techniques employing interferometers 
(Reference 1) concluded that they could not be used, due to the extreme phase dis- 
tortions caused by multipath reflections. These distortions prevented correct 
ambiguity resolution of a long baseline interferometer with a short baseline inter- 
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ferometer, since the short baseline unit did not have enough accuracy to determine 
the correct grating lobe of the long baseline unit. 

By using an array of interferometer pairs, and recursively resolving ambiguities 
from a shorter baseline to a slightly longer baseline, correct ambiguity resolu- 
tion can be obtained for arrays of any length. Further reduction in multipath 
caused errors can then be obtained by averaging the results of Interferometers 
at different positions. The MADGE system is currently configured as a Category II 
landing aide and, thus, is not designed to meet all of the ILM requirements. This 
analysis examines the feasibility of extending the design concept to meet Category 
III requirements . 

Summary of Capabilities 

Space diversity has been shown to be an effective technique for the reduction of 
specular multipath errors in an interferometer. The MADGE elevation array shows 
errors of less than 10 mrad even in severe multipath conditions when operated as a 
focused array. The current MADGE system is apparently operated in an unfocused 
manner, but range estimates to allow focusing could be made available with minor 
change in system operation. The 10 mrad error exceeds the allowable 1.2 mrad error 
for an MLS elevation unit, but is close enough that proper site selection should 
yield an acceptable system. 

Two fixes were attempted. 

o Stuffing the array with additional antennas, keeping the 
spacing between antennas of a pair constant. 

o Extending the array by 4 more antennas, maintaining the array 
spacing constant and thereby increasing the spacing between 
antennas of a pair by a ratio of 3/2. 

Stuffing the array did not yield any improvement in performance. Extending the 
array provided a slight Increase in accuracy, bounding the error below 7.5 mrad 
rather than 10 mrad. Thus, neither of these particular fixes are sufficient to 
provide Category III performance in a severe multipath environment. 
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The MADGE azimuth unit does not utilize space diversity for multipath protection. 

In the cases examined, which were severe multipath situations, the azimuth unit 
had errors large enough to cause ambiguity resolution errors. The azimuth unit 
does have some multipath protection, and would probably be adequate in a less 
severe environment. It would be better to use a space diversity array similar 
to the elevation unit to insure operation in any multipath environment. 

In the box canyon environment, an elevation type array operated as an azimuth 
unit showed better performance, holding peak errors to 13 mrad rather than the 
18 mrad of the normal MADGE azimuth unit. More importantly, the space diversity 
array is capable of providing correct ambiguity resolution in the presence of 

multipath distortions of this magnitude, whereas the standard MADGE azimuth array 
is not. 

In all cases using a space diversity array, large multipath errors tend to be highly 
oscillatory with a mean error much smaller than the peaks of the oscillation. 

Since no filtering is employed in the computation of these multipath errors, it is 
likely that the errors could be reduced to a much more acceptable value if filter- 
ing were employed. For example, in the case of the extended elevation array at a 
1° glide slope in the Samoan airport environment, the peak errors are about +7.5 
mrad, whereas the mean error is .02 mrad. This is certainly an acceptable multi- 
path error component. 

The MADGE equipment as it now exists is not capable of Category III operation. Severe 
multipath environments can cause errors to exceed the limits for CAT III equipment, 
focusing errors cause excessive geometric biases, and the system is not configured 
to provide roll out guidance. The following changes should be made to MADGE to 
provide a workable Category III equipment; 

o The interferometer arrays should be focused . 

o The azimuth array should be replaced with a space 
diversity array. 

SYNTHETIC CIRCULAR APERTURE SYSTEM 
Concept Description 

An independent landing monitor concept called ASTROLABE is being developed by ONERA 
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in France. In this system the two angular coordinates, elevation and azimuth angles, 

. I 

are obtained by 2-dimensional phase modulation of a microwave carrier, induced by 
the continuously rotating motion of a pair of antennas. 

A conventional interferometer array consisting of two elements separated by several 

f 

wavelengths has high angular resolution within an ambiguous lobe structure created 

-:fi 5 

by the grating lobe pattern. Resolving the ambiguity consists of determining which 
of the many possible lobes contains the real target. Rotating the array provides 
ambiguity resolving capability since the grating lobes are not evenly spaced in 
angle. The cross-correlation of the received signal parameters with the theoretical 
signal parameters based on the rotating geometrical grating lobe structure provides 
a high-precision angle determining capability. The scheme is conceptually simple 
but difficult to apply because of the complexity of the 2-dimension cross-correla- 
tion processing requirement which is characteristic of synthetic aperture systems. 
Optical processing techniques are used in many synthetic aperture systems because 
of the 2-dimensional transform characteristics of lens systems. Although digital 
techniques are applicable, system costs are much higher because of the fast 
sampling speed, high computation rates, and large mass storage requirements of a 
near real-time system. ONERA personnel have patented a rather simple and compara- 
tively economical optical system that provides a high precision signal processing 
capability that is applicable for real-time display of angle proportional information 
for either ground or aircraft installation. Current displays used in the develop- 
ment tests give a visual display of the location of the correlation point whose 
position is representative of the angular coordinates referenced to the rotating 
array . 

Several different system concepts using the Astrolabe principle are possible and 
are discussed in Section III. Two concepts were considered representative and 
examined in more detail in the analysis: variation 1, an Astrolabe rotating 

receiving array on the aircraft with several beacon transmitters highlighting 
the runway extremities to provide a real-time perspective display of the runway 
extent and variation 2, an Astrolabe rotating receiving array on the ground with an 
aircraft beacon and a two-way DME/information data link providing near real-time 
cockpit indications of aircraft approach corridor azimuth and elevation angular 
coordinates and range-to-go. The two variations are shown in Figures 3 and 4. 
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4 BEACON TRANSMITTERS 
AT RUNWAY EXTREMITIES 
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Summary of Capabilities 

The concepts used in the Astrolabe system show potential for use as an ILM. 

Reference 1 contains the desired ILM angular accuracy specification of 0.10° to 

0.12° for 20 accuracy. To assure 2cr accuracy of 0.12° with Astrolabe would require 

element separations of approximately R/X=80 which would require spinning element 

separations (R) of approximately 3 meters at X-band and 0.8 meters at K^-band. The 

selection of frequency for the transmissions is somewhat arbitrary although it 

should be recalled from Reference 1 that approximately (1.3 x R (Km)) dB more 

max 

power is required for K -band operation than X-band operation in heavy rain (16mm/ 

sl 

hr). For example, guidance or position determination at a range of 8 Km would 
require over 10 dB (10 times) more power at K^-band than at X-band for the same 
receiver signal to noise ratio. 

It is conceivable, however, since an 80 cm element separation would be possible 
in the nose of a commercial aircraft, that K -band operation of a spinning antenna 
in the aircraft for high accuracy angle determination would be feasible. This 
configuration would have the requirement for stabilization of the spinning antenna 
to eliminate high frequency roll, pitch, and yaw angular deviations to achieve the 
required high accuracies. A gimballed system for the spinning antenna (with the 
spinning elements generating the gyroscopic stability platform) could be used with 
sensors to derive smoothed aircraft to antenna reference angles to relate derived 
angles to the aircraft reference frame. 

Multipath errors consist of two types: direct multipath terms which are unresolv- 

able from true targets over short observation times and cross-product multipath 
terms that originate from the product detector which are resolvable over relatively 
short (less than one second at X-band) integration times. Analysis shows the cross- 
product terms are reducible by at least 12 dB but the direct multipath terms are 
left to the expertise of the display analyst for removal. The multipath problem 
therefore favors ground receivers and analysts since the multipath geometry would 
be more familiar and predictable. A pilot would see a somewhat different display 
to interpret at each airport. 

Comments on ILM applications up to this point have been general for nearly any 
configuration. When multiple systems are used in an effort to attain location by 
triangulation or geometrical shape orientation, several other considerations are 
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necessary. 


Multiple ground receiving locations were shown to require large separations (2 to 
4 Km) . These wide separations seriously hinder displaced ground receiving sites 
since the triangulation function must be performed centrally. Either multiple 
receivers, receiver processors, and analysts or multiple receivers with IF or 
detected IF signals relayed to a central processing and analysis center would be 
required resulting in a costly configuration. 

Configurations with airborne receivers and displaced ground transmitters have the 
same problem of wide separation of ground antennas for high accuracy although the 
analysis is centralized to the aircraft. The wide separations necessary for ranging 
accuracies necessarily limit the achievable angular accuracies. This effect 
coupled with multipath problems will degrade performance such that any multiple 
configuration will be of questionable use for flare, touch-down, and roll-out 
information. 

The multiple ground antenna configurations, although capable of providing high 
accuracy results prior to flare maneuvers are not practical ILM's because of the 
complexity of the analysis, interpretation, and triangulation requirements. If 
the high accuracy requirements can be waived it would be feasible to use variation 
1 (Figure 3) where the imagery (4 beacon points) would be representative of the four 
corners of the runway and presented to the aircrew in the form of a pseudo real- 
world perspective image for a VFR type of approach. This application would require 
the spinning (or a stepped Wollenweber array) receiving antenna in the aircraft. 
Decreased accuracy requirements would negate requirements for a stabilized antenna 
platform. Beacon location (point) accuracies on the display to the order of 1° 
should be achievable. Suitability of a four-point image of the runway for guidance 
to the pilot was not determined in this study. 

A high accuracy ILM application at a reduced cost when compared to the multiple 
ground antenna configurations would consist of a combination of a single beacon 
transmitter on the aircraft, a spinning receiving antenna on the ground, and an 
interrogator/response TOA (time-of -arrival) ranging capability incorporated into 
the ground- to-alrcraft data link for relaying angular information to the aircraft. 
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This system would provide high precision approach corridor information (azimuth 
and elevation angle along with precise range-to-go as described in reference 1, 

ILM Multilateration Sensor section). It should be pointed out that this combi- 
nation of sensors would benefit from the most accurate capabilities of both 
techniques: Astrolabe precision angle referencing without triangulation 

techniques and its associated accuracy degradation and high precision ranging 
by TOA multilateration techniques without GDOP (Geometrical Dilution of Precision) 
inherent in multilateration angle determination. This concept should be capable 
of providing 2 sigma angular accuracies of 0.1° and ranging accuracies to a 
few meters which are consistent with ILM requirements. 


CONCLUSIONS 

1. The use of a large thinned array receiving antenna in the airborne bistatic 
radar concept provides a significant resolution improvement over a comparable 
monostatic radar. Part of this improvement can be traded off for the ILM 
problem through reduction of operating frequency from Ka-band to X-band with 

a corresponding drop in atmospheric signal attenuation. This allows the system 
to operate at reduced output power and, thus, reduced transmitter costs. This 
cost advantage is offset, however, by the cost of the wing-mounted phased array. 
Total system cost of the X-band bistatic radar and the Kg-band imaging radar 
would be •approximately the same. Since runway detection range of both radars 
is approximately the same due to similar beamwidths, the advantages of the X- 
band bistatic system become subjective; quality of function (weather radar and 
ILM) and the availability of more equipment expertise since X-band technology 
is more developed than K^-band technology. 

These advantages are at the expense of increased installation design costs. 

From preliminary examination, it appears feasible to install the thinned array 
antennas in the wings of most commercial aircraft. However, each aircraft type 
will probably require its own unique installation design. 

2. The complex interferometer/DME combination (MADGE) shows a potential of meeting 
Category III accuracy requirements for ILM, however, as currently configured 
(designed for Category II operations) it does not meet the accuracy requirements. 
Analysis of the elevation antenna array and signal processing showed this concept 
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can suppress multipath Interference and provide angular accuracies equivalent 
to MLS to within 200 to 300 meters aircraft to ground antenna ranges. Use of 
the elevation concept for both azimuth and elevation angle measurements and 
addition of a second set of antennas to provide flare, touchdown and rollout 
guidance should satisfactorily meet the Basic and Case I ILM requirements. 

Use of radar altitude beyond runway threshold for flare and location of the 
azimuth antenna at the stop end of the runway should eliminate the need for 
the second set of antennas . 

3. Analysis of the circular synthetic aperture concept (Astrolabe) showed azimuth- 
elevation angular measurements of sufficient accuracy (0.1°) could be achieved 
over the desired ILM operating range. Use of multiple beacon configurations to 
measure range through triangulation were not satisfactory and it was concluded 
DME must be Included in the concept to provide an independent range-to-gc 
measurement . 

Electro-optical processing is used to cross-correlate the received signal with 
a reference pattern to generate a pseudo-real world display of the azimuth- 
elevation direction to the signal sources from the receiver. Multipath sources, 
interference rings and other noise sources also appear on the display with 
human interpretation required to identify the desired signal source. We con- 
cluded this type of display may be feasible for a ground controller where each 
aircraft approach would generate a similar looking display. A display of this 
type in the aircraft would show a different noise pattern at each airport 
requiring longer response time by the operator for interpretation. Further 
signal processing may also be feasible to electronically measure azimuth/eleva- 
tion angles with output to a more common guidance display. 

4. Comparative costs for the three systems in this report and the feasible ILM 
concepts of Reference 1 were generated internally through comparisons to like 
equipments. They are based on 1974/1975 costs and do not represent specific 
manufacturer inputs. These costs and weights of the airborne units are as 
follows : 
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Airborne Cos t 


Grd. Cost 


Airborne Wt. 


Bistatic Radar $62K - 82 kgms 

Complex Interferometer/DME $13K $41K 8 kgms 

Circular Synthetic Aperture $15K $23K 8 kgms 

(variation 2) 

RECOMMENDATIONS 

1. We can not recommend a specific sensing concept which stands above all others 
from our analysis in meeting the ILM requirements. We have identified several 
concepts, including the three in this report and six in Reference 1, which show 
potential capability to meet at least part of our baseline ILM requirements. 

Each system is costly and/or will require major development efforts before 
becoming operational in addition to each having unique problem areas ranging 
from the human interface to installation difficulties both on the ground and 

in the aircraft. 

Radars, operating in the X to Kg-band, have potential for meeting to some degree 
each of the major ILM functions examined; approach, landing and rollout guidance, 
obstacle detection, and real world imagery. With reduced environmental extremes 
and possible enhancement of runway/background differences, development of use- 
ful operational systems should be feasible. 

Several concepts; multilateration, redundant MLS, complex interferometry, 
synthetic circular aperture, provide greater approach coverage and higher 
accuracy for the ILM functions of approach, landing, and rollout guidance than 
the radars but provide no capability for obstacle detection or real world 
Imagery. Selection of one of these concepts as an ILM requires reliance on 
another set of equipment or operational procedures to cover the obstacle avoid- 
ance function and verification of crew acceptance of the non-real world guidance 
display (cross-pointers, synthetic image, etc.) for the confidence building 
function. 

2. As discussed in the recommendations of Reference 1, there is no clear consensus 
of opinion in the industry on the functional usage or purposes of an ILM. This 
study started with the baseline ILM function of landing monitor for a Configura- 
tion K jMLS -equipped aircraft. Using this function, the associated ILM sensing 
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requirements are generated by the possible Category III environmental and 
operational landing conditions. 

We were unable in this studj; due to time allotted, to explore in detail other 
possible functions of an ILM or frequency of occurrence considerations of the 
various environmental /operational conditions leading to extreme sensing require- 
ments. These factors must be analyzed and traded off with the potential 
sensing concepts before an optimal concept can be selected for any user. 
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SECTION III 
DETAIL ANALYSES 


BISTATIC RADAR APPLICATION ANALYSIS FOR ILM 
Introduction 

A bistatic radar uses separate transmitting and receiving antennas such that scat- 
tered incident energy may be received at angles other than the Incidence angle , 
Receiving antennas incapable of handling large transmitter powers may be used, or 
arrays may be thinned with resulting grating lobes being removed by selection of the 
second antenna. This application is of the latter type, using a special phased 
array receiving system (Reference 2) . This application has several advantages 
which will be presented later. 

The previous analysis of Independent Landing Monitors (ILM's) using radar 
techniques (Reference 1) concluded that X-hand radiation had the potential for 
adequate penetration of Category III weather conditions (very short or "zero" 
visibility) but did not have the capability for high resolution because of the 
large physical size of a monostatic radar antenna. This bistatic application 
allows a large physical receiving antenna for high resolution at X-band by locating 
array elements (individual antennas) along the lower edge of a wing. The receiving 
function is simplified by "thinning" the array (removing array elements) such that 
fewer array elements are used while maintaining the long array length for high 
resolution. Fewer array elements (16 to 32 elements) make it possible to perform 
the receiver heterodyning at the array elements and thus minimize transmission 
line losses and simplify routingof the array signals through the wing structure 
by replacing the RF waveguides with IF cabling. Thinning the array does create 
the problem of ambiguous grating lobes of the receiving array, however, by properly 
selecting array element spacings to match the transmitting antenna, it is possible 
to place the receiving grating lobes in the nulls of the transmitting radiation 
pattern. This concept is the application analyzed by this report. 

A representative scenario for the ILM application is an aircraft during straight- 
in approach to a 4000 meter runway during adverse weather conditions (Reference 1) . 
The aircraft is on a 3° glideslope with a 5° nose up attitude relative to the glide- 
slope and a velocity of 75 m/sec (145 kts). The analysis describes the expected 
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signal power ratios as a function of : 

Range 

Weather Conditions 

Runway and Surrounding Culture Scattering Cross Sections 
Equipment Parameters 

The specific parameters that were used in the analysis for each of the above factors 
are described in the following paragraphs. 

Analys is 

The analysis of the bistatic radar application was performed with the same consider- 
ations that were used in the previous ILM analysis so that comparisons could be 
made with previous results. The same analytic software routines were used for 
deriving single-pulse return power ratios and the additional factors considered in 
Reference 1 were reviewed. The single-pulse results and the additional factors 
are discussed in detail in Paragraph C below. The parameters and environmental 
factors considered in the software analysis are defined and reviewed in Paragraphs 
A and B respectively. 

A. Radar Parameters 

Two specific radars were analyzed as illuminators for a bistatlc application. 
Parameters typical of a Marconi LN-66 radar were used because of previous 
bistatic tests of this very short pulse, high PRF radar (Reference 2). For 
comparison, a stock RCA AVQ-30X weather radar as an illuminator was analyzed 
since this radar is already available on some aircraft. The parameters used 
in the analysis are summarized in Table 1. Table 2 lists the receiving system 
parameters. The array and receivers as described in Reference 2 were analyzed 
due to their apparent applicability for aircraft mounting. 

B. Software Scenario 

The software previously developed (see Reference 1 for description) was used 
for the single-pulse return power calculations. The basic four most adverse 
weather cases of each type of weather (2 fog, 1 rain, and 1 snow) were used to 
determine the radar penetrability. The four cases are summarized in Table 3. 

The X-band backs cat taring cross-sections, in dBSM (dB referenced to one 
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Table 1. Typical Radar Transmitting Parameters 


(References 1 and 2) 


Radar System 

Marconi LN-66 
(BISTAR System) 

RCA AVQ-30X 

Band 

X-band 

X-band 

Frequency 

9.4 GHz 

9.4 GHz 

Power 

12 kw 

65 kw 

PRF 

5000 

200 

Pulsewidth 

40 nanoseconds 

5 microseconds 

Antenna Gain 

33 dB 

35 dB 

1st Sidelobe (H) 

-21 dB 

-30 dB 

(rel. to main lobe) 



1st Sidelobe (V) 

-21 dB 

-30 dB 

(rel. to main lobe) 



Az Beamwidth (one-way) 

2.3° 

2.9° 

El Beamwidth (one-way) 

3.6° 

5.2° 

Scan Parameters 



Scan Sector 

o^ 

o 

o 

1+ 

o 

o 

180" (±90°) 

Scan Rate 

I200°/second 

45°/second 

Images Per Second 

20 images/second 

1/4 image/seconds 
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Table 2. Receiving Array Parameters 
(Reference 2) 

Array Element 

Flared Slot 

Max Gain - 20 dB* 

Az Beamwidth 60° (see Figure A-2a) 
El Beamwidth 7° (see Figure A-2b) 

Array, linear 

Max Gain* - 5.8 dB + 10 ^ ^°®10 

= 20.2 dB 

( j = 7.5, so K' - 6.5) 

Number of Elements, n = 32 
. Array Length 7.7 m (25 feet) 

Az Beamwidth 0.23° 

El Beamwidth 7.0° 

1st Sidelobe (H) , -30 dB 
No sidelobes on the Elevation Pattern 


1 

Transmitter Pulsewldth (sec) 

*See Appendix A for Derivation 


Receiver 

Noise Figure 9 dB 
Noise Bandwidth, B, Hz = 
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Table 3. Analytical Weather Cases 


Weather Case 

Description 

Density or Fall Rate 

Altitude 

D 1.3 

Heavy Rain 

16 nm/hr 

ground to 1500 m 

D 2.4 

Radiation Fog 

4g/m3, zero visibility 

ground to 60 m 

D 3.4 

Evaporation Fog 

o 

4g/m , zero visibility 

ground to 1500 m 

D 4.3 

Heavy Snow 

9 mm/hr 

(equivalent water) 

ground to 3000 m 


Table 4. 


X-Band Backscattering Cross-Section, 



dBSM 


Asphalt 

-56 

to -45 

Concrete 

-61 

to -38 

Grass 

-42 

to -5 

Snow 

-43 

to -29 


- 29 - 












square meter) are listed in Table 4. The extreme values for grass backgrounds 
and a 23 meter wide runway were used to derive the results in Paragraphs D.2. 
The remaining parameters and expressions are identical to the previous 
analysis . 

C. Analytical Considerations 

The computer runs provided the single-pulse results shown in Tables 5 and 6. 

The very low backscattering cross-sections assumed for the results of Table 5 
prompted the use of circular polarization for precipitation clutter reduction 
(25 dB reduction for rain and 16 dB reduction for snow). Additionally a modi- 
fication of the Marconi LN-66 was run to test the feasibility of degrading 
resolution to increase grass background returns. This modification, designated 
XRADAR, has identical parameters to the LN-66 except the pulsewidth was 
increased (decreasing the receiver bandwidth) by a factor of 5 (to 200 nano- 
seconds) and the receiving array was shortened to 16 elements (12.5 foot length, 
beamwidth of 0.46°). This resulted in a decrease in array gain of 6 dB, an 
Increase in footprint and clutter volume by 10 dB and a decrease in receiver 
noise bandwidth by 7 dB (1/5). 

Several additional factors were considered as modifiers of the single-pulse 
results to account for statistical and dynamic properties of complex multiple- 
pulse returns and the display or monitoring functions. The factors considered 
in the previous ILM sensor analysis were reviewed and the terms appropriate or 
significant to the bistatic scenario are discussed in this section. 

% 

Lj - Scanning Loss 

The losses due to scanning of the antennas are of three origins. 

Lj-|^ - loss incurred because the scan rate is too fast or 
the PRF is too low to receive at least one pulse 
return from each resolution element. 

Lj 2 ~ loss due to mis-tracking between the two separately 
scanned antennas. 
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Table 5. Radar Return Power Ratios for a smooth concrete (y = ~61 dBSM) runway surrounded by 
short dry grass or snow (y = -42 dBSM) (circular polarization for precipitation 
clutter reduction is assumed) 
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Table 6. Radar Return Power Ratios in dB for a rough concrete runway (y = ~38 dBSM) 
surrounded by 8" to 10" vegetation (y = -5 dBSM) 
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Lj 3 “ loss from electronically scanned phased array, the gain 
decreases as approximately |cos 0| for 6 measured from 
broadside. 

For a beamwidth of 0.23°, the Marconi radar as presented in Reference 2 has been 
modified to transmit one pulse per resolution element and 20 independent pulses 
for each resolution element (one per scan) per second allowing the possibility 
of direct CRT display. (Independence will be shown later.) The RCA radar has 
a sufficiently slow scan that one pulse is radiated each 0.225° or approximately 
once each resolution element of the 0.23° array. This is coincidentally con- 
venient, however, only one return from each resolution element occurs every 4 
seconds. This data is too slow for direct viewing without display storage 
such as a scan converter or storage tube in the display. Is thus negligible 

although the effects of the slow sample rate will be evident in other factors. 

Mis-tracking between the two antennas is not a serious factor since both 
illuminating (transmitting) antenna beamwidths are at least a factor of 10 
wider than the phased array. Reference 2 and Appendix A demonstrate that 
weighted amplitudes and/or non-linear element spacings can reduce sidelobe 
(and grating lobe) levels such that spurious returns can be suppressed to at 
least 30 dB below the receiving array mainlobe. As long as the phased array 
can track the illuminating antenna within approximately 1°, Lj2 will be neg- 
ligible. 


The loss due to electrically scanning the phased array 0 degrees from broadside 
is approximately |cos 9|. At the scan limits, this would be about a 1.5 dB 
reduction in gain. Large crab angles necessitated by high cross-winds may 
require operation in this region. Other effects that would require consider- 
ation in this region would be degraded beamwidth (beamwidth is proportional 

to -1 — Q -r ) and array length across which the propagating short time-duration 

I cos y I 

return pulse must propagate. The receiving array length and transmitter pulse- 
width for this analysis were acceptable. (e.g. a 40 nanosecond pulse could be 
detected by nearly all elements of a 25 meter length array at 6 as large as 


30° ( 


3 X 10® m/sec x 4x10“® sec ~ 


sin 30 


= 24m). L ]^3 is negligible for targets 
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directly ahead and approximately a 1.5 dB loss for targets at the scan limits 
of the Marconi radar. Operation of the RCA system beyond 30° would cause 
larger losses, however, operation at such extreme angles is not envisioned 
for ILM's at this time. 

L_ - Sensor Losses due to aircraft motion and instabilities. 

a 

The high resolution (0.23°) capabilities of the receiving array will dictate 
that stabilization is required since rms angular yaw rates of approximately 
0.8 degrees/sec with deviations of several degrees were Identified in the 
previous study (Reference 1) . Tracking of the yaw instabilities within 
approximately the beamwidth (0.25°) is possible with on-board instrument sensors. 
Pitch and roll sensor inputs are similarly required. L is thus assumed negli- 
gible^ provided display stabilization is incorporated using aircraft instrument 
sensors . 

Lj - Scan Distribution Loss 

The derivation of scan distribution loss will follow the same assumptions as 
the previous analysis; specifically, the bistatic ILM display will occupy a 
small amount of aircrew workload time resulting in only a couple of seconds 
for aircrew viewing during each required reference of the ILM display. The 
display integration is thus limited to a few seconds and the number of pulses 
(samples) and independent scans forming the reference display image are 
bounded by those few seconds. The Marconi illuminator may provide approxi- 
mately 40 samples of each resolution element in 40 scans (20 scans per second) 
while the RCA system would provide only one sample (if at all) in the same 
time interval. In view of this rather severe restriction, we'll assume the i 

scan limits of the RCA radar are modified to +30° at the same scan rate (a 
relatively minor modification) such that up to two samples may be received 
in the same two second time interval. Recall from the previous study that scan 
distribution loss is derived (from Reference 3) by holding the cumulative 
probability of detection, P^,, constant in the expression: 

= 1 - (1-Pd)^ (1) 

where P^, - cumulative probability of detection in j scans 
Pjj - single scan probability of detection 
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The variation of single scan probability of detection, Pj, as a function of 
increasing number of scans, j, results in a decreasing S/N requirement. As 
the number of scans Increases, the number of pulses available for integration 
in each scan decreases, offsetting the decreasing S/N required. The effect of 
the trade-off between single scan and multiple scans is shown in Figure 5 . 

The scan distribution loss L^j is relative to direct integration of all pos- 
sible pulses (assumed correlated within one scan and uncorrelated between 
successive scans) falling within the Interval of consideration. 

Lj - Fluctuating Signal Loss 

The fluctuations of the pulse returns contribute losses, Lj, as a function of 
number of pulses Integrated and probability of detection. Swerllng and 
Schwartz have characterized the effects and the loss as a function of proba- 
bility of detection is shown in Figure 6. 

Lci - Collapsing Loss 

Four types of collapsing loss are considered in the bistatic system: noise, 

Lcn; beam, Ld2» I'ange, and display, 

The noise collapsing loss, LcH, arises from noise sources within parallel 
receiver channels adding coherently. Sources would Include effects such as 
local oscillator noise at the mixers of parallel receiver channels to the 
phased array adding at IF. The effects are reducible by circuit techniques 
(e.g., balanced mixers) and minimizing the coherency of the noise prior to 
summation by adjusting the noise path lengths as described in Reference 2. 
Careful adherence to proper noise limiting techniques should allow reduction 
of L^ii to negligible levels. 

The beam collapsing loss, Lj,22» is considered in part by the analytic software 
beam considerations which generate beam patterns representative of t3rpical 
patterns to add in contributions from sidelobes and grating lobes. Other 
effects such as beam pattern distortions during operational flights due to 
effects of flaps or landing gear deployment can only be speculatively consid- 
ered pending operational tests of the actual array. Assume the array is mounted 
as far out on the wing as possible to remove the nose wheel and engine nacelles 
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as far as possible from the array scan coverage and as far back on the wing as 
possible to minimize pattern blockage by theleading edge flaps. This config- 
uration should minimize beam collapsing loss, T-‘d 2 > resulting in only a slight 
increase in precipitation clutter return caused by "knife-edge" diffraction 
along the lower edge of the flaps. This diffraction distortion of the antenna 
pattern would tend to increase the vertical beamwidth (collecting more precip- 
itation clutter due. to increased vertical dimension of the clutter cell), 
however the nearly vertical distortion would not degrade resolution since the 
vertical angle reference is established through ranging and knowledge of air- 
craft altitude and glide slope. Since the precipitation clutter will be shown 
to be substantially less than background returns, the effect of this loss 
should be small. 

The range collapsing losses, returns from ambiguous ranges at high 

PRF is expected to be severe due to the extreme variations in backscattering 
cross-sections. The previous analysis suggested variable PRI’s as a solution 
for integrating out this effect. Reference 2 also suggests pulse-to-pulse 
beam agility although this implementation would require an electronically- 
scanned phased array illuminator or several transmitting antennas. Variable 
or staggered PRI is recommended to minimize to acceptable levels. 

the display collapsing loss is minimized by providing appropriate cockpit 
Instrument sensor inputs so that effects of yaw and roll image distortions can 
be minimized to the appropriate sensor accuracy. Typical gyro/servo accuracies 
of nearly 0.25° would be in close agreement with array resolutions. 

L. - Integration Loss 

The integration loss (relative to n integration gain) is a function of the 
integration efficiency of the display/display processing and the pulse return 
correlation functions. Figures 7 and 8 from Reference 1 provide the correla- 
tion information for background and precipitation returns. The display is 
assumed to integrate correlated pulses with an efficiency of n^*^ = integration 
gain (Reference 1), uncorrelated as n^'^. The correlation data shows that the 
background and precipitation clutter are essentially uncorrelated from pulse 
to pulse as the standard deviation of doppler velocity within 10 km is within 
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Figure 8. T^, Time Between Samples 



1.5 m/sec which is the region of inseparable combined turbulence, shear, and 
ground wind doppler velocity components. 

D. Results 

The overall results of the single pulse and multiple pulse factors define a set 
of complex system requirements needed to minimize the loss factors presented 
previously. A summary of the loss terms is presented in Table 7 and Figure 9. 
The multiple pulse effects are insignificant compared to the magnitude of the 
individual pulse power ratios. Tables 5 and 6 are summarized in the graphs 
in Figures 10 through 13. The graphs have curves labeled to show the most 

significant mechanisms contributing to the curve. Thus, if Bg>>S and 
C<<N the curve approaches Bg/N. The different slopes evident in the curves 
are usually indicative of differing dominant terms at near and far ranges. 

The AB reference curve is derived from the convoluted power return maximum 
variation from the runway and background as the runway is scanned. Thus, 
in order to "detect" the runway the range must be less than the range at 
which the AB reference curve crosses the +6dB minimum detection threshold. 

The analysis showed the application has potential for providing high resolution 
direct imaging of the runway during approach in severe rain, snow, or zero 
visibility fog. Specific results are Itemized below: 

1. Lowest backscattering cross-sections provide inadequate contrast, 

AB, and background to precipitation clutter, Bg/C, to detect the 
runway unless circular polarization is used for clutter reduction 
and degraded resolution is used to decrease the noise bandwidth. 

With these assumptions, approximate detection rangesxfor 20 dB 
contrast (runway to background) are: (XRADAR results) 


Weather 

Detection Range, 23m Runway 

16 mm/hr Rain 

(1 . 9 km) 

9 mm/hr Snow 

(2 km ) 

fog, zero visibility 

(1.9 km) 
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Table 7. Summary of Multiple Loss Factors 
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Figure 10. Summary of Power Detection Ratio 

for Marconi LN-66 Bistatic C-4^ 
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Figure 11, Summary of Power Detection Ratio for 
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Marconi Bistatic System (Maximum Y ) 
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Figure 13. Summary of Power Detection Ratio ^ 
AVQ-30X Bistatic System (Maximum Y) 


2. Runway images derived by passive reflector enhancement or imaging 
of cultural objects such as runway lights or other objects with ap- 
proximately 1 square meter scattering cross-sections can be detec- 
ted with the extended detection ranges shown: 


Weather 

Detection Range, 23m Runway 

16 mm/hr Rain 

(2.6 km) 

9 mm/hr Snow 

(3.0 km) 

Fog, zero visibility 

(2.8 km) 


I 


3. Higher backseat tering cross-sections allow further utilization of the 
high resolution capabilities of the Marconi illuminator with the 
array receiving system. For this case, a runway to background con- 
trast of 35 dB was analyzed. The approximate detection ranges are: 


Weather 

Detection Range, 23m runway 

16 mm/hr Rain 

(4.1 km) 

9 mm /hr Snow 

(5.4 km) 

Fog, zero visibility 

(4.4 km) 


4. The limitations above for the 23 meter runway are due to the runway 
to background contrast. This contrast of 35 dB is approximately 
the maximum allowable based on the 33 dB sidelobes measured for the 
array (Reference 2). Using this as the limiting case, assuming 6 
dB Bg/C minimum required for detection, the maximum precipitation 
rates may be derived for the detection of a 23 meter runway at 1 n.mi. 
and 2 n.ml. These maximum rates are derived from the approximate 
expression C a R^*^, where R is the precipitation fall rate in mm/ 
hr. (Reference 1) 
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Maximum Precipitation Rates for 23 meter Runway Detection 


1 NM 

Rain 

67 mm/hr 


Snow 

59 mm/hr Equivalent Water 

2 NM 

Rain 

20 mm/hr 

Snow 

17 mm/hr Equivalent Water 


COMMENTS ON ILM APPLICATION 

The application of the bistatic concept as a ILM would require a system 
consisting of a short pulsewidth, high PRF illuminator similar to the 
Marconi LN-66 with a phased array receiving antenna with narrow beamwidth 
resolution, essentially the system proposed in Reference 2. Due to high 
resolution (small resolution cell) capability of the system, extensive dis- 
play processing (such as with a digital scan converter) would not be bene- 
ficial since the complexity imposed by the large number of possible image 
locations would make a digital application very time consuming and costly 
in mass storage. The correlation functions of the signal returns do not 
indicate the potential for significant power ratio enhancement by utilizing 
such a processing capability beyond that which can be achieved with direct 
CRT display and nominal phosphor persistence to minimize flicker. 

The bistatic system is shown in block diagram form in Figure 14. 

The physical equipment requirements of this system are large requiring modi- 
fication of existing radar equipment (if a suitable radar is available as 
an illuminator) or installation of a suitable radar illuminator. Installation 
of a receiver, installation of a cockpit display subsystem, and installation 
of a large external receiving array with internal routing of signal infor- 
mation. The system obviously requires an extensive aircraft retrofit. 
Additionally, since the array is rigidly mounted to the aircraft, a "tie-in" 
to motion sensors is required to stabilize the high precision imagery. The 
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Figure 14. Bistatic Radar ILM Application 
















use of existing instrument sensors or the addition of instrument motion 
sensors would be necessary for accuracies of less than 1°. The list of 
equipment/modifications necessary for high accuracy utilization of a BISTAR- 
type system are listed below with their functions. 


Transmitter (illuminator) (i.e., modified LN-66) 

Narrow Scan - high sample rate 

Narrow Pulsewidth - high resolution 

High PRF - short range, high sample rate 

Staggered PRF or Pulse-to-Pulse - interference suppression 
Frequency Agility 

Narrow Beamwidth (a few degrees)- interference and sidelobe suppression 


Weather Radar Modification (i.e., AVQ-30X modification) 


Decrease Scan Limits 

Decrease Pulsewidth 

Increase PRF (coordinated 
with decreased pulsewidth to 
provide the same duty cycle) 

Staggered PRF or Frequency - 

Agile Transmitter 

Receiver 

L.O. Locked to Transmitter 
Heterodyning performed at Array - 
Fast Receiver Blanking - 

Manual and/or Automatic Gain, 
Threshold, and Linearity Controls 


increases sample rate 
decrease clutter cell size 
increase sample rate 

enhances interference suppression 


maximize received signal 
minimize RF loss 

minimize nearest effective range 
and protection of components 

optimize receiver to scenario 


Array, Receiving (at least 20 to 30 elements) 

Electronically Phase Scanned - highest scan rates and resolution 
and Focused 

Wing Mounted (underside) - maximize array length 
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Directional Elements 


- higher gain and interference 
immunity 

Thinned Array Element Spacing - suppresses receiving pattern 

matched to Transmitting Antenna grating lobe ambiguities 

Display and Processor 

Direct CRT Display for High - near real-time 

Pulse Rates 

Perspective Display - real-world similarity 

Stabilized Display - more accurate viewing, less "jitter" 

and "ghosting" 

Mutual interference effects are reduced by using staggered pulses or a 
frequency agile transmitter and the highly directional antennas. The very 
small backscattering ciross-sections (-40 to -60 dBSM) typical of highly 
specular surfaces (e.g., concrete, asphalt, snow) will necessitate a large 
dynamic range for the receiver such that direct interference signals do not 
disrupt the reception of following cultural returns. With sufficient dyna- 
mic range, the cultural returns will integrate to effectively reduce the 
sporadic interference effects. 

The production costs of the bistatic radar (Figure 1) are based on use of a 
standard weather radar, such as the RCA AVQ-30X with small modifications for 
adding a second mode of operation, and the addition of the thinned array 
antenna. Projected cost of the X-band radar with modifications is $30,000. 
Cost of the thinned array antenna is $1,000 per element including phase 
shifters and phase shifting logic. For the 32-element antenna, this 
corresponds to $32,000. Total cost per aircraft for this ILM concept is 
estimated at $62,000.00. This concept is totally contained on the aircraft 
with no additional ground system costs. Estimated weight of the airborne 
system is 82 kgms . It should be noted that approximately $25,000 of this 
figure could be attributed to the weather radar mode. Aircraft installation 
costs were not estimated. A major cost add-on for this system is the non- 
recurring installation design costs for each specific aircraft type to 
install a thinned array antenna in the wings. 
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ANALYSIS OF SPACE DIVERSITY INTERFEROMETERS FOR ILM 


INTRODUCTION 


A space diversity interferometer utilizes a number of interferometer pairs 
spaced different distances from the multipath sources. Since multipath- 
caused errors in an interferometer can be positive or negative depending on 
the path length difference between the direct and multipath signal, averaging 
the results over multiple interferometer pairs can yield an increase in 
accuracy over that obtainable with a single Interferometer. 

The particular application being analyzed here is the MADGE system, developed 
by the MEL Equipment Co. Ltd, Sussex, England. It utilizes a space diversity 
interferometer to obtain elevation data, a standard Interferometer with a 
special array ambiguity resolution technique to obtain azimuth, and a pulse 
DME to obtain range. The previous analysis of Independent Landing Monitor 
techniques employing interferometers (Reference 1) concluded that they could 
not be used, due to the extreme phase distortions caused by multipath 
reflections. These distortions prevented correct ambiguity resolution of a 
long baseline interferometer with a short baseline interferometer, since 
the short baseline unit did not have enough accuracy to determine the correct 
grating lobe of the long baseline unit. 

By using an array of Interferometer pairs, and recursively resolving 
ambiguities from a shorter baseline to a slightly longer baseline, correct 
ambiguity resolution can be obtained for arrays of any length. Further 
reduction in multipath caused errors can then be obtained by averaging the 
results of interferometers at different positions . 

Since the previous analysis showed that specular multipath was the major 
problem with interferometry, it is the only area addressed in this analysis. 

Various flight paths and terrain scenarios were used to establish different 
multipath situations. Specific parameters for the various situations are 
given in following discussion. 
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ANALYSIS 


The primary effort in this analysis was directed toward tile creating of a 
computer program general enough to handle the MADGE space diversity array 
and to create multipath situations which could challenge its capabilities . 
The previous analysis (Reference 1) showed that only specular multipath 
prevented interferometry from being effective, so specular multipath was 
the only factor considered in this report. 


A detailed description of the MADGE concept is contained in Reference 4. 

Only key factors will be repeated here. 

A. MADGE Parameters 

The MADGE system consists of three interferometer arrays, one for azimuth 
and two for elevation. 

The azimuth array consists of seven horn antennas, spaced such that a 
pair of antennas is available with spacings increasing in a 2:1 ratio 
from wavelength to 16 wavelengths. 

The azimuth angle data is derived completely from the 16 wavelength 
long interferometer. All the other antenna pairs are used to resolve 
ambiguities. Assuming a high signal to noise ratio, ambiguities can 
be correctly resolved assuming multipath phase distortion remains less 
than 60° with a 2 to 1 ratio interferometer array. 

The main elevation array consists of eight antennas arranged in a linear 
vertical array with 3.6 wavelengths between antennas. Since the main 
array has a 5 degree grating lobe width, a smaller array is also provided. 
It is a linear array of four antennas with a 1.8 wavelength spacing 
providing unambiguous angle measurement to 25.6 degrees elevation. 
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INDEPENDENT ANGLE ESTIMATES 



OVERALL ANGLE ESTIMATE 


Figure 15. Effective Elevation Array Connection 


The elevation angle data is the average of the estimates from 4 pairs 
of interferometers as shown in Figure 15. Figure 15 does not represent 
the actual connection of the equipment, but rather the manner in which 
the information is used. The array is divided into 4 pairs, each pair 
with a separation of 14.4 wavelengths between antennas. Ambiguity 
resolution is accomplished first by using the data from the antennas 
as a recursive array using the 3.6 separation to resolve the 7.2, then 
using the 7.2 to resolve the 10.8, etc., and by using an additional 
array with smaller spacing to resolve the ambiguities in the 3.6 
wavelength interferometer. This technique will correctly resolve 
ambiguities in the presence of multipath phase distortions up to 90°. 
Assuming this condition to be met, each of the four pairs provides an 
independent measurement with different multipath. The average of the 
4 measurements may therefore be better than any individual measurement. 


Other equipment parameters are given in Table 8. 


B. Software Scenario 

The software described in Appendix B was exercised for fifteen different 
situations (Table 9) . The first twelve are related to the elevation 
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TABLE 8 


MADGE SYSTEM PARAMETERS 


Band 

C-Band 

Frequency 

12 channels from 5-5.25 GHz interferometer 
for angle measurement, pulse DME for range 
measurement, ground derived and data linked. 

Power 

200 watts 

Bandwidth 

3 MHz 

Minimum Operating SNR 

15 dB 

Size 

Aircraft Unit 


2 - 3/8 ATR short 
8 Kg 


Ground Unit 


Azimuth Antenna 


1.5 M X .3 M X .4 M 
38 Kg 


Elevation Antenna 


1.9 M X .4 M X .3 M 
40 Kg 


Transponder 


•6Mx .3Mx .4M 
23 Kg 
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Table - 9 Scenario Test Situations 
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Figure 16. Flat Surface 




Figure 17. Samoa Airport 
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array and the last three are related to the azimuth array. Each run is 
divided into two segments, the first from a range of 5000m to 1000m, the 
second segment from 1000m to 20m. 

Runs 1, 3, 4, and 5 were done using a rough flat surface, of permittivity 
80 + j 40 and surface RMS deviation of .05 wavelengths. Run 2 had no 
multipath and is a check on residual errors. Runs 6 through 12 use the 
actual approach path to Pago Pago airport in Samoa, and is an attempt to 
provide an extreme multipath situation. Run 13 is a hangar scenario 
where the aircraft approaches between rows of hangars, and is the situation 
which MEL claims to have tested its equipment in. Run 14 is an imaginary 
box canyon or fjord, with two aircraft on a taxiway in front of the system. 
Both run 13 and run 14 are extreme multipath situations. Run 15 is the 
box canyon environment again, using a space diversity array in place of 
the MADGE azimuth unit. 

Three pages of output are provided for each run. The first page provides 
the setup of the system and the environment. The data titled "ANTENNA 
PAIRS" is the height of each antenna above the reference plane in meters. 
FLIGHT PATH SEGMENTS provides the starting and ending coordinates of the 
aircraft in meters relative to the reference coordinates. On those pages 
where the setup for two runs is made, the first two columns of flight 
path data refer to the first run number, and the second two columns refer 
to the second run number. Any surfaces which could cause multipath have 
their parameters listed under the heading "REFLECTING SURFACES". The 
second page of the output is a plot of angular error in radians versus 
range in meters from 5000m to 1000m. The third page is a similar plot 
from 1000m to 20m. These plots are automatically scaled by the computer 
routine, and thus the long range and short range data usually have 
different scales. 

Analytical Considerations 

The program developed for this analysis computes the effect of specular 



reflections on a space diversity interferometer of up to ten pairs of 
antennas, for up to ten flight path segments, with up to ten reflecting 
surfaces. The program description is given in Appendix B. 

Simple large target reflection is all that is considered. At each 
increment of each flight path segment for each surface and each antenna, 
a computation is made to determine whether or not a point exists on the 
surface such that the angle of incidence from the aircraft equals the 
angle of reflection to the antenna. If such a point exists, the rough 
surface modified complex Fresnel reflection coefficient is computed, 
and used in the determination of that surface's contribution to multipath 
errors. If no such point exists with the bounds of the surface, no 
multipath contribution from this surface (for this antenna at this 
increment on the flight path) is computed. All surfaces are assumed to 
lie normal to the plane of interferometer operation, i.e., for the 
azimuth interferometer, all surfaces are assumed vertically oriented and 
of infinite extent in the vertical direction. In the plane of operation, 
surfaces are represented by a line segment. 

Diffraction is ignored, except that no test of blockage of a surface 
reflection by another surface is performed. Hence, some diffraction 
effects are included unintentionally. 

The analysis mechanized in the program is basically the technique des- 
cribed in Reference 1 under the section Azimuth Channel Multipath, with 
a few modifications. Logic is included to average the errors over an 
array, and actual reflection coefficients are computed. Since a range 
estimate is available, provision for focusing of the interferometer 
array is included as an option. Focusing compensates for the spherical 
nature of arriving wavefronts, and the displacement of the individual 
interferometer pairs from the reference datum and each other. 

Reflected signal strengths are modified by the ratio of the direct range 
to the antenna from the aircraft to the range via the reflecting surface. 
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the magnitude of the reflection coefficient, and the ratio of the antenna 
pattern in the direct direction to the pattern in the direction of the 
reflected ray. A very simple approximation to the antenna pattern was 
used, since the only available data was the 3dB beamwidth. The patterns 
used are shown in Figures 20 and 21. Polarization is assumed to be 
vertical relative to the surface of the earth at the point of reflection. 

A total vector sum of the field strengths at each antenna is computed, 
and an estimate of physical angle is made, using a focused or unfocused 
algorithm as specified by the data and assuming ambiguity resolution can 
be properly obtained. This angle is compared to the actual (physical) 
angle, and the resulting error is output. 

D. Evaluation (See pages 71 through HI for computer runs) 

Run 1 

The first run is the basic MADGE elevation array measuring an aircraft 
on a 3° glide slope in an environment consisting of a single flat re- 
flecting surface. The results show that multipath distortions are 
smoothed quite well, but errors of up to 7.5 mrad are still present. 

These errors occur at short ranges of 200 m or less. Thus, with careful 
siting this could be acceptable. 

Run 2 

The focused algorithm was run without any multipath distortions to 
determine the residual errors inherent in using range to the antenna 
rather than range to the reference point and errors due to truncation 
in the computer. The results show errors on the order of a few micro- 
radians, so the program was judged acceptable. 

Run 3, 4 and 5 

The next three runs are the focused MADGE elevation array operating in 
a flat reflecting surface environment, measuring aircraft on 1°, 3°,, and 
6° glide slopes for runs 3,- 4 and 5 respectively. 
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Figure 20. Antenna Field Pattern Plot 
Vertical Array 
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Figure 21 . Antenna Field Pattern Plot Horizontal Array 



The results are as could be expected in that at very low elevation 
angles which correspond to high reflection coefficients and inability 
to obtain any reduction in multipath signal strength via antenna patterns, 
the errors are highest, running about 10 mrad at 1°, 2.5 mrad at 3°, and 
less than 1 mrad at 6°, all at long range. 

The runs at 1° and 6° show an increasing error at short range. These 
runs were offset 55m to the side of the elevation array however, and 
the magnitude of the error is similar to the difference between the 
elevation angle from the reference datum at the base of the interferometer 
array and the elevation angle from a reference point at the intersection 
of the ground with the flight path. That this is a purely geometrical 
effect is shown by the 3° run, which has no offset. In the 3° run, the 
errors decrease with decreasing range as the incidence angle to the 
reflecting surface becomes greater and the antenna pattern becomes more 
effective in reducing the multipath signal strength. 

Runs 6, 7 and 8 

The runs use the MADGE elevation array in the Samoan airport scenario at 
elevation angles of 1°, 3°, and 6°. Although a 1° glide slope is not 
possible at this airport, those surfaces which would interfere with a 1° 
glide slope are excluded from reflecting by the specular angle logic 
included in the program, so that for the 1° case the situation is 
effectively an airport with a 10m cliff at the approach end of the runway. 
In all of these runs, the aircraft flight path is offset 55m from the 
elevation array horizontally. 

In run 6, the long range data shows a small error. The multipath 
indicators show that only surface number 2 is reflecting, so this is to 
be expected. In the short range data, both surface 1 and surface 2 are 
reflecting and oscillations in the multipath error of ±10 mrad are 
observed in the data. 
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Runs 7 and 8 show oscillations in the multipath errors on the order of 
1.2 mrad, which is still within the specification. In addition, the 
geometrical effects of the offset are again observed. 

Runs 9 and 10 

Runs 9 and 10 are' again the Samoa airport at 3° and 6° glide slopes. The 
interferometer used here is a MADGE array, stuffed with 8 more antennas 
(4 more pairs) in between the antennas of the basic MADGE. The results 
are almost Identical to runs 7 and 8, so this is not an effective 
modification. 

Runs 11 and 12 

For runs 11 and 12, four additional antennas were added to the top of the 
elevation array, and antenna pairs are defined as every sixth antenna 
rather than every fourth antenna as in the existing MADGE system. Run 
11 is a 1° glide slope and Run 12 is a 3° glide slope. Some improvement 
in accuracy can be noted relative to runs 6 and 7. 

Run 13 

Run 13 was made with a single interferometer pair in azimuth. This is 
effectively the configuration used by the MADGE system, as all the 
remaining antennas are used only for ambiguity resolution. The multipath 
scenario consisted of 4 hangars near the interferometer and a truck. 

At long range, the data is very good on the order of microradians. The 
multipath indicators show that none of the objects are influencing the 
measurement. At short ranges, various objects are reflecting over part 
of the path, and oscillations of ±18 mrad are shown in the data. However, 
these oscillations correspond to phase distortions of ±100 degrees . 

Since, for proper resolution of ambiguities, the MADGE azimuth array 
cannot tolerate phase errors greater than ±60 degrees, the MADGE azimuth 
unit could not operate in this scenario. 
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Run 14 


Run 14 was also made with a single interferometer pair, with a scenario 
representing a box canyon or fjord. 

The results are similar to run 11, except that the multipath problem 
exists at all ranges . Rapid oscillations in measurement error of ±18 
mrad occur over the entire flight path. As with run 11, errors of this 
magnitude cannot be accommodated by the ambiguity resolution technique 
so that actual error would be much larger, on the order of several 
degrees due to Improper ambiguity resolution. 

Run 15 

Run 15 is an Identical environment to run 14, but the basic MADGE elevation 
array is turned on its side to use as an azimuth unit in place of the 
regular interferometer MADGE uses in azimuth. 

Oscillatory errors on the order of ±13 mrad exist in the data. The 
significant difference from Run 14 is that the space diversity array 
can function properly with errors this large, whereas the regular inter- 
ferometer cannot. 

RESULTS 

Space diversity has been shown to be an effective technique for the reduction 
of specular multipath errors in an interferometer. The MADGE elevation array 
shows errors of less than 10 mrad even in severe multipath conditions when 
operated as a focused array. The current MADGE system is apparently operated 
in an unfocused manner, but range estimates to allow focusing could be made 
available with minor change in system operation. The 10 mrad error exceeds 
the allowable 1.2 mrad error for a MLS elevation unit, b^it is close enough 
that further effort on the array would probably yield an acceptable system. 

Two fixes were attempted. 

o Stuffing the array with additional antennas, keeping the spacing 
between antennas of a pair constant. 
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o Extending the array by 4 more antennas, maintaining the array 
spacing constant and thereby increasing the spacing between 
antennas of a pair by a ratio of 3/2. 

Stuffing the array did not yield any improvement in performance. Extending 
the array provided a slight increase in accuracy, bounding the error below 
7.5 mrad rather than 10 mrad. Thus, neither of these particular fixes are 
sufficient to provide Catetory III performance in a severe multipath environment. 

The MADGE azimuth unit does not utilize space diversity for multipath 
protection. In the cases examined, which were severe multipath situations, 
the azimuth unit had errors large enough to cause ambiguity resolution 
errors. The azimuth unit does have some multipath protection, and would 
probably be adequate in a less severe environment. It would be better to 
use a space diversity array similar to the elevation unit to insure operation 
in any multipath environment. 

In the box canyon environment, an elevation type array operated as an azimuth 
showed better perfoinnance, holding peak errors to 13 mrad rather than the 
18 mrad of the normal MADGE azimuth unit. More importantly, the space 
diversity array is capable of providing correct ambiguity resolution in the 
presence of multipath distortions of this magnitude, whereas the standard 
MADGE azimuth is not. 

In all cases using a space diversity array, large multipath errors tend to 
be highly oscillatory with a mean error much smaller than the peaks of the 
oscillation. Since no filtering is employed in the computation of these 
multipath errors, it is likely that the errors could be reduced to a much 
more acceptable value if filtering were employed. For example, in the case 
of the extended elevation array at a 1° glide slope in the Samoan airport 
environment, the peak errors are about ± 7.5 mrad, whereas the mean error 
is .02 mrad. Errors on the order of .02 mrad are certainly acceptable 
multipath error components. 
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SNUMR = Gina'S. ACTIVITY f = 01, REPORT CODE = 12, RECORP COUfoT = 000270 
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RUN 1 
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PLOT OF ANGULAR ERROR VERSUS RANGE 
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RUN 2 



PLOT ANGULAR ERROR VERSUS RANGE 
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RUN 5 Segments 3 & 4 
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PLOT OF ANGULAR ERROR , VERSUS RANGE 


m 

O 



I o 

I 

I Ui 

I o 
I o 

I ^ 

I • 


I 

I 

I o 

I 

I UJ 

I o 

^ I o 

I 



I 

I 

I fo 

I o 

I 

I u. 


CO 

u 

0) 

4J 

M 


fg 

O 

LU 

C9 

oo 

UJ 



o 

< 

ID 

o 

q: 

cc 


I rg 
» O 
I 

I UJ 

I c 
— < I o 

I c 

I • 

I o 

I 

I «M 

I o 

I 

I UJ 

I o 

^ I o 

I o 

I • 

I 

I 

I «NJ 

I o 

I 

I Uj 

I o 
— I o 
I o 
1 • 

I <v 


; X /V rg m 

0 o o 

1 I I 

UJ U.I UJ 

o o o 

o o o 

m o o 

• • • 

^ 


en 

O 

I 

UJ 

o 

o 

o 


rg 

O 

I 

U; 

O 

O 

o 


rg 

O 

I 

u. 

o 

o 

(TV 


rg rg fg 
C O c. 
I I I 
UJ U1 L. 
O O O 
O O O 
O m O 


rg 

C 

I 

U.' 

o 

o 

tn 


o 


(suBipey) yoyya 


79 



SNUMB = G1005« activity <f = Cl* RFPC’T COPF = I?, PfCCR'' CCUM = 000270 


O 


o o o o o o 


o o o o o o 


o o o o o o 


*Ti 

U.I 


n 

o 

or 


I 

iij 

o 


o 

o 


o o o o o o 


O O t"5 O O O 


oooocooooo 


>- o 

> •••••••••• 

— oooooooooo 


or 

Ijj 

a 


o 

u.' 

O 


O 

o 


c a c cj o c 


ocoooooooo 


o o o c- o o 


«/} 

ui 


• • 




o 

o o 



LJ 





O' 





■ o 




o 

o 

>» 




u 

c 



• 

O’ V 

u 

•£> ^ 

0 

0 

0 

0 

D 

o 


o 

r- o 




o 

• • 



c 

o 

t>« *»4 

m f\j rsi 

t— t 

z 

c 


o O O 

o 

u. 



U: U- a. 

L. 


o 


o mo 

tf’ o 




O' rj O 

O IP 




<7* u> r«* 




or < 


• 



in a* 

0 

0 

0 

0 

0 

o 



• • 





^ O 

ir. 


4T 



t~ mm 

CM • 

Ul NJ 



<: c CO 

O 

h- 



Ul 


< 



2r. Uj tU 

b 

^ ■ 



O- r\J O 

m o 




Ll O o 

f\j ip 

c 

• 

o CD 

m rvj <7 

IP 

O' 

c 

^ *o 

• • • • • 

• 

c 


• « 

X o o o o o 

o 

c 


^ o 


lO 

iP o 



< 

»- 

UJ 

m 


u. 

z 

u o 

o 



Ui 

< z 



*- 

z 

u_ •« 

u 


z 

UJ 

cr 

o 


ir. 

ft 

r> o: X 

o 

<M o 


\J 

ip <■ 

o 

(M m 


z 

1- 

o 

« • 

U XX 

I ^ 

C 1/) 


^ O 

H- h- Wl tO 

%n 

z 

• - 


K- a a. cr — 

•— u. 


O ' 


U. < < < X z 

z o 

»- UJ 

1 




u u a 


a; or 

Q:«/)t/^cOLL.u-U.a: 

u; < Ui 


UJ UJ 


UJ 

^ U- (D ^ • 

Q. 2 


►-.cr 

u. o: z 


a o 


UJO o 


3 U 

< 

17 

or IP Z 



X > M X >- fSl Z 




•— oooooooooo 


COCCOOCrCOO 


oooooooooo 


- 80 - 


RUN 4 



PLOT OF ANGULAR ERROR VERSUS RANGE 



OOOE 02 1.500F 03 2.500F 03 3.500F.03 <i.500E 03 

l.OOOE 03 2.000E 03 3.000E 03 <>,OOOE 03. 5.000E 03 



PLOT OF ANGULAR ERROR VERSUS RANGE 



(sue^ppy) yoyy3 


RANGE (Meters) 
RUN 4 



PLOT OF ANGULAR ERROR VERSUS RANGE 


I 

( 


I 

I 






tr. 

o 

I 

U- 

o 

c 

o 

in 


m 

c. 

t 

L.- 

O 

o 

lA 

• 

rg 


O 


lA 


(suBipey ) yoyy3 


83 . 


l.OOrF 0? ?.000E 03 3.000F 03 <>,000E 03 5.000E 03 

OCOE 02 1.300E 03 2.500E 03 3.500E 03 4.500E 03 



PLOT OF ANGULAR ERROR VERSUS RANGE 


m 

o 


(suB^pey) yoyy3 


RANGE (Meters) 
RUN 5 



o 


o o o o o o 


o o o o o o 


o o o e o o 


M 


X 

i9 


O O O O 
till 
Ui UJ UJ Ui 

o o o o 
o o o o 
o o o o 
o o o o 
M «> «r\ 

oooooooooo 


o o o o o o 


o o o o 

UJ teJ 111 Ui 

o o o o 
o o o o 
o o o o 
o o o o 

^ ^ 


«« oooooooooo 


o 

o 

o 


o 

s 

\J 

UJ 


o 

o 

\J 


• • 

o o 


o o o o o o 


o e o o o o 


o o o o o o 


« 

UI 

& 


M ru 
o o o o 

UJ UI Ui UI 

o o o o 
o o o o 
o o o o 
o o o o 

«P o o o 


oooooooooo 


(N m 

o o o 

Ui Ui UJ 

o o o 
.000 
000 


o 

s 

u 


*0 

>000000000 


CM 


eO 


1—H 

CO 

t/) 

c/> 

4-> 


C 

c 

0) 

0) 

£ 

E 

o> 

o> 

Q) 

Q> 

00 

00 


00 

z 

z 

ZD 

ID 

ai 

0 ^ 


^ 

0000 



0 ^ 

0 

0 

0 

0 

0 

0 

X 

0000 


>• 

0 







0000 


ft 

• • 






0 

0000 


0 


m 





X 

0 « ^ 


ft 

UJ 

ft 


0 


0 0 


0 

UJ 

ft m 



UJ 


Ui Ui 


Ui 


oooooooooo 



ft 


m 0 


ft 0 




« 


ft 


<u 0 


0 ft 






ft 


ir\ t >4 






0 

• 'C 

• 

• 

• • 

• 

• 


CM r 4 «M 




0 

0 

000 

0 


000 


n 

• • 

•4 0 

«A 





ft 

UI UJ Ui 




K ft 


m m 


fU 

UJ isj 

000 




ftO 


00 


0 

h- 

000 


>- 


UJ * 





ft 

000 




ft UI 


UJ UJ 


UJ 

X 

0 ft r- 




0 ft 


ft 


m 0 

M 

^ 0 ft 


> 


UJ ft 


ft ft 


N ft 

0 


• • 


0 • 

U) ft 


«u ft 


ft 

ft 

oooooooooo 



• 

• 

• • 

• 

• 

0 

1 


u 

• • 

xo 

0 

00 

0 

0 

0 



< 

«>«o 

*- 




ft 

ft 4 J 





< 




P~ 

Ui 

fM ft ft 


• 


a 




X 

u 0 

000 


lU 






Ui 

ft X 



e 

ft 

*- 




X 

u. 

Ui UJ UJ 


0 


X 




UI 

p:»- 

000 


UJ 

< 

49 




ft 

5 ft X 

000 


0 

a tu 0 





u 

ft ft 

000 



IM »\ 





X 

K 

0 «0 


M 


U. 


X 

X 

X — 

49 ft 

ft rg ft 



ft-«o 

K 

¥- 

^ «n 

«o 

ft 

X 


• • 

CD 

ft 

►- ft 

ft 

ft 

•m 


M 

dodboobo 

0 0 

S 

Ui 

luft 

ft 

ft X 

X 

X 0 

h* Ui 



5 



►- 




U U ft 



ft 

ft ft ft 

ft 

«A 

u. 

u. 

U. ft 

Ui ft UJ 




<UiUj 

U 




Ui 


ft 0 


a B 

ft K 


^ H- 


K CD 

U. ft X 




a 0 

•«< 

ft 

ft ft 

ft 

ft X 

UI 3 9 




0 -1 

ft 




3 

ft ft X 





X 

>• 

N X 


N X 





- 85 - 



PLOT OF ANGULAR ERROR VERSUS RANGE 



RANGE (Meters) 
RUN 6 



PLOT OF ANGULAR ERROR VERSUS RANGE 


o 



(suBipey) yoyya 

- 87 - 


RANGE (Meters) 

run 6 



SVUM0 = DFD02, ACTIVITY f = 01, REPORT CODE = 12, RECORD COUNT = 000270 
ANTENNA PAIRS 


o 


O O o o o o 


o 


o o o o o o 


I/) 




in 

^ ,-4 ^ 


o 

(U 

o o o o 

• • 


2 

1 1 1 1 

o o 


Z 

tli UJ liJ UJ 


o o o o d o 

19 

O O O O 


n 

O O o o 



o 

o o o o 



a 

o o o o 

rv «A (A tA 




oooooooooo 

• • 

o 



o o 

o o o o o d 


<V rg <NJ «M 

O O O O 




UJ UJ Ul UJ 




o o o o 
o o o o 



>- 

o o o o 




o o o o 





o o 

o 

> 





oooooooooo 



M 



o o o o o d 





a 

04 €SI C4 



Ui 

o o o o 



a 

ui Ul ui UJ 

• • 



o o o o 

o o 

o 


o o o o 
o o o o 
o o o o 

00 CQ 0> CO 


o o d d d d 


oooooooooo 


<v CM m 
o o o 


ix a 
UJ UJ 
Q. -k 

a o 

3 -J 


^ yf 
O 


CD yC 
ITt <D 


O CD 


o 


UI UJ UJ 

o o o 


«A 

o o o 


UJ 

o o o 

o d d d d d 

H- 

lA r- m 


< 

>0 O' ^ 


2 



O 

QC 

O 

oooooooooo 

o 

o 

CM ^ 


u 

o o o o 


19 X 

Ul UJ UJ UJ 

o d o o d d 

2 

o o o o 



o o o o 


c 

o o o o 

m Aj <M ^ 

2 

O CD ^ ^ 

o o o o 

UJ 

m CM m m 

UJ UJ Ui UJ 

C7* ro O lA O 

O' <M O O VN 

O' JA rg ymt 


ooooododod 

• ••••• 


CM AJ <V 

o o O o o o 


o o o 

) 

lA 

UJ UJ Ul 

• <A EA <M 

Ui IS 

o o o 

: o o o o 


O O o 

1 

< 

o o o 

: Ul UJ UJ UJ 

2 

o *A r- 

•O' <M O «A o 

M 

'O O' 

1 O' ^ O' 4M lA 

o 


AJ O' lA 

a 

oooooooood 

• ••••• 

o 

1 

: o o o o o o 

o 


4A 

40 U 



UJ 

<M <4- -4^ 

2 

U 19 

o o o 

UJ 

< 2 


X 

u 

Uj UJ UJ 

UJ 

a *- 

o o o 

O' 

O O' X 

o o o 

SJ 

«A < 

o o o 

2 

K 

O CD ^ 

z z z ^ 

19 4A 

fo CM m 

»- K >- tA V) aA 

2 


a a a ^ ^ u 


oooooooooo 

<<<2220 

UJ 


i-i 

u u oc 


iA iO tA u u u. or 

UJ < UJ 


Ui 

_JU<D^fNm^JA <0 CD o o 

CD 

u cr z 


<<<<<< 2 

Ul o => 


O 

X V KI K >> rsl 2 

<rio 2 



88 


RUN 7 



PLOT OF ANGULAR ERROR VERSUS RANGE 



(suBipey) yoay3 


89 . 



PLOT OF ANGULAR ERROR VERSUS RANGE 



(sueipey) yoay3 


90 


OOOE 02 3.000E 02 5.000F 02 7.000E 02 9.00QE 02 

2.000E 02 4, OOOE 02 6.000E 02 8.000E 02 I 1*000E 03 



PLOT OF ANGULAR ERROR VERSUS RANGE 



RANGE (Meters) 
RUN 8 



PLOT OF ANGULAR ERROR VERSUS RANGE 


m 

o 





o 


O O O 0-0 o 




o 






6 o o o o o 

lA 






lA 

m 4 Mrf 




o 

bJ 

o o o o 



• • 


z 

• 1 1 1 



o o 


X 

lb lb lb lb 





19 

o o o o 




o o o o o o 

3 

o o o o 





O 

o o o o 





QC 

o o o o 
rg m m m 






oboQoooooo 



• • 

o 





o o 

6 6 o 6 b o 


eg eg eg eg 
O O O O 






bJ UJ lb lb 
O O O O 
O O O O 



rg O 


>■ 

O O O O 



00 mU 


K 

O O O O 



• • 


MM 




mH m4 

o 

> 






•M 

oooooooooo 













o 



d d d d d o 

X 






QC 

eg eg eg eg 





b/ 

o o o o 

oO 

ir 

'<1’ <M 


Q. 


o 

O 9‘ 



lb lb lb lb 


o 

• • 



o o o o 

o» 

o 

*-i o 

o 


o o o o 
o o o o 

CO 

II 




o o o o 






2H 

K 





3 

z 


o o o o o o 


oooooooooo 

QC 

3 





o 

•43 <4h 





u 


m (M m oj <b ^ 





• • 

o o o o o o 


eg eg m 

CO 

c. 

O 



o o o 

or 


111 LlJ ll.i U1 IaI lAJ 



s: 

o 


tf-. o 1 ^ a o ^ o 


bJ b b 

o 

u 


o> in o ^ in o 1 ^ 


o o o 


lU 


O' in MM mm (M 

lA 

o o o 

1— 

a 


• ••••• 

bJ 

o o o 



c o o o o o 

►- 

in e- 

t—t 

• 

OD -O 


< 

«0 O MM 

Q 

rg 

^J iTi 


z 




• • 

^ rg m K> eg fO 

M 

oboooobboo 

o 

o 

II 

M- O 

o o c o o o 

a 

a: 




Ui Ui Ui Ui bi Ui 

o 



Ui 


e> o m o m o 

o 

eg <# ^ 


o 


c m rg o- tn O m 

u 

o o o o 


c: 


•4- m m m M^ 




u 


• • • • • . • 

e X 

b‘ b‘ b- b 




o o o o o o 

z 

o o o o 



r*- o 


MM 

o o c o 


<y 

• • 


o 

o o o o . 



Ml4 

e' eg isj eg eg mm 

r 

O m ^ ^ 


a 


O O O O O C 

bJ 

CA m lA 


IX! 






0: 


Lj lx; bJ bj b; b^ 

O O ro o o in o 


oooooooooo 




O' m eg o in o in 

O' m m eg m M* 




o 

cn <• 

• ••••• 


f\i ej eg 



(TN CD 

o o o o o o 


o o o 


II 

• • 

MM O 

lA 

lA 

b b b 


% 


K nj lA <e. rg eg 

UJ eg 

o o o 




2000000 


o o o 


>- 


b: 

< 

o c o 


► 


S' b,' b.' b: bJ lx b-j 

z 

one- 


MM 


OO OegO 0 1^0 

MM 

mM o 


> 


b o IT' 'C o m eg IT 

c. 



MM 

O OD 

lA ^ in eg o in in 

a 

oooc oooooo 


K 


• ••••• 

o 

1 



• • 

X o o o o o o 

c- 



< • 

MM O 

K lA 

lA b 





< h- 

lb 

<M 


* 


Cl 7. 

U V9 

o c o 



»/> 

b> 

< z 



C 

ct 

1- >. 

b ^ 

b b b 


C> 

MM 

X bJ 

a 1- 

o o o 


MM 

«■ 

i.': o 

D O X 

o o o 


O 

Q rg o 

o 

b'. < 

o o o 



«n 

w z 

h~ 

o CD mm 


II 

< • • 

b I I X — 

O lA 

fC* eg lA 



i: ^ o 

i/i «/•. i/> 

•2 





»-cratt»-MM*-bL 

MM 

OC OOOOCOOO 


> 

b. 

b<«.<X'Z20 

►- b- 




h- 

MM MM l-M 

U u o: 



?• 

cr fr 

o•4/:«/■^«nb.bbnr 

b 4 bi 




< Uj UJ 

u ui 

^ b a 

— •rg<A'*rin^T''a c o 



a 2 

fr H- O'. 

b o. 7- 

mM 



a o 


b X X 




:d -j 

< ~ 

CY tj 7’ 





X >- fv X V k; 7 


-93- 




PLOT OF ANGULAR ERROR VERSUS RANGE 


O 



(suBipcy) yoay3 


94 



PLOT OF ANGULAR ERROR VERSUS RANGE 


fn 

o 



(sueipey) yoyy3 


95 


RANGE (Meters) 
RUN 9 



PLOT OF ANGULAR ERROR VERSUS RANGE 



l.OOOf: 03 ?.000E 03 3.000F 03 ^.OOOE 03 S.OOOE 03 

5.000E 02 1.500E 03 2.500E 03 3.500E 03 ♦.SOOE 03 5.500E 03 



PLOT OF ANGULAR ERROR VERSUS RANGE 




o 


o o o o o o 


o o o o o o 


O O O O O o 


0 o o o 

1 I I I 
Uj ui llJ UJ 

o o o o 
o o o o 
o o o o 
o o o o 

fM tfN 

oooooooooo 


o o o o o o 


(N rsj ro rg 
O O O O 

UJ UJ UJ Ui 

o o o o 
o o o o 
>■ o o o o 

^ o o o o 

> ««•••••••• 

M 00.0 0000000 


O O O O O O 


o o o o o o 


m M 

o o o o 

Uj W 

o mo o 
O* <M o o 

9> M 

• ••••• 

o o o coo 


ft <Ni fg 

UJ o o o o 

a 

iU OJ UJ UJ 

o o o o 
o o O O' 
o o o o 
o o o o 

00 <0 03 03 

oooooooooo 


CVJ ^ 

oS oQ 

T— I CO 

in m 
-M -M 
c c 
0 ) O) 
E E 

CJ> CD 

a> a; 

00 CO 


# 

O rg 




«C ft 

rg 

m rg 



z 



• • 


m tn fg 


oooooooooo 


rg ^ 

O 

O o o 

o 


II 




ft 




UJ 

UJ Um Ui 

o 


Ui 


ft 

rg ft mo 

c 

rg ft ft ft 

c 


ft 

*0 ft rg u> 


O O c o 

c 


ft 

fsj ft 



u 


• 

• • • • • 

C X 

U U b Ui 


CM '•t 

o 

o o o o o 

ft- 

O O O O 

»- 

— o 




O O C O 

ft 

• • 



c 

O O O O 

c 

rg 

ft 

rg rg o 

ft 

O « — ft 

a 


o 

O O O 

u 

m rg m m 

Ui 






ft 


U' 

U^ u. u- 


O C C O O O C O o c 



O 

me ft c 



♦ 


o 

r- o ft m 





m 4 

^ rg m 



o 


• 

• • • • • 


rg rg rg 


ft 00 

o 

o o o o c 


O O O 

II 

• • 






^•o 

iT 


1 /; 

Ui UJ U 


T o o o o o o 


ooococooco 

i 




o 

►- 















ft 







»- 

UI 

rg 


ft 

tr. 



CL 







ft- 

UI o 

o 

& 

o 










u: 

ft. :? 




o 



K- 







3 

u 

u; 

U 

u> 




I 







UJ 

ft ►“ 

o 

c 

o 




K.0 







ri 

“ ft- 

>. o 

c> 

C: 

a 

oc 

O 

« 







U 

4 ft ^ 

o 

c 

O 


ir 

m 








Z 


o 

CO 


<. 

• 

• 

u 




X 

X 

X 


o ^n 

m 

rg 

fT» 



c 


»- 

K- 




%/ 


z 

• • 

t 

• 

ft 



»- 

a 

a 

o 




u. 


0 

0 

O 

o 

u; 



u 

ft 

ft 

ft 

Z' 

ft 

if 

c 

»- u, 




►- 



ft 

h- 

H- 

X- 

M 




U W ft 




ft 

ft 

ft 

sf; 


(/• 

u 

u. u. 

ft 

Uj ft 

UJ 



< 

iLf UJ 

u 







UJ 

^ u 

ff: ^ rg 

m 

ft 


a. 

ft 

K- 

►- 

I— 

»- 

►- 


ffi 

u. rv 

3 l 




a 

o 


< 

ft 

< ft 

ft 

ft 

?: 

Ul=) 

D 




98 


PLOT OF ANGULAR ERROR VERSUS RANGE 


m 

o 



RANGE (Meters) 
RUN 11 



PLOT OF ANGULAR ERROR VERSUS RANGE 



OOOE 02 6.000B 02 l.OOOF 03 1.400E 03 1.800E 03 

4.000E 02 8.000E 02 1.200E 03 1.600E 03 2.000E 03 



PLOT OF ANGULAR ERROR VERSUS RANGE 



(suBippy) yoay3 


101 


OOOE 02 1.500E 03 2.500F 03 3.500E 03 ^.500E 03 

l.OOOE 03 2. OOOE 03 3. OOOE 03 4.000E 03 5. OOOE 03 



PLOT OF ANGULAR ERROR VERSUS RANGE 



(SUB^pey) yoyya 
- 102 - 



AIRCRAFT FLIGHT PATH 


.o 


o o o o o o 


-o o o o o o 


o o o o o o 


eg 

0 o o o o 

1 I I I I 

t‘j iij 4 U 

O o o o o 
o o o o o 
o o o o o 
o o o o o 


oooooooooo 


o o o o o o 


^ ^ 

o o o o o 

lU UJ U1 ID Ul 

o o o o o 
o o o o o 
o o o o o 
o o o o o 


oooooooooo 


o o o o o o 


eg fM rg eg eg 

o o o o o 


2 

D 

o 

o 


o o o o o o 


o o o o o 
o o o o o 
o O O 0.0 
0.0 o o o 

eg eg eg eg eg 

oooooooooo 


o 

u 


o o o o o o 


Id ID ID ID ID 

o o o o o 
4/1 o o o o o 
ID . O o O O o 
*- lr^ O 'O o o 

< ^ ej eg 

z •••«•«•••* 

oooooooooo 


o 

o 

SJ 


o o o o o o 


^ eg eg eg eg o 

o o o o o o 

ID ID Id ID Id Id 

o o irt o o o 
o ir> e> o le' A# 4 TW 

^ lo eg in «ei 

• ••••• 

o o o o o o 


2 . ID (D ID ID ID ID 

ooof^oomo 
Lu o r* ooNin 
4/1 lO eg OB ^ 

• ••••• 

o o o o o o 


XIX*- 

4 /) 40 40 

o: a or *-*-•« It 
< < 4 z z z o 

^ ^ ^ M Mi M 

4/1 4/) 40 u. Id u. tt 


<<<<<< X 
D 

K >- eg X >• fsi z 


O 
o 
u 

O X ID ID ID U; ID 

~ o o o o O 

O O O O O 

O W> O O O 

» 4 O 

eg ^ it ^ 


Z — 4 


OOOOOOOOOO 


) ID Uj UJ ID ID 
J eg o o O O O 
O O O O O 

: o o o o o 

4 / 40^00 
« — « it eg m 4 

> •••••••«•• 

oooooooooo 


o 
o 

in SJ 

ID eg eg ei e> w 

4D O O O O O C 

< z 

ID — lD< U< Ui Id ID 

CC K O O O O C 

noixocooo 
t/i < o o o o o 

t- O 4T* it — O 

0 4/1 o* UT m it m 

z ••»••••••• 

«« OOOOOOOOOO 
4- ID 
U U tt 
ID < ID 

_ilD(D'-fgeiiti/'^e~«c-oo 
Id <r X ^ 

ID3 D 
O' 4/1 Z 


103 


RUN 13 



PLOT OF ANGULAR ERROR VERSUS RANGE 


O 



(sueipey) ypyi)3 


104 . 



PLOT OF ANGULAR ERROR VERSUS RANGE 


I 

I 


• I 

I I 




I 

I 

I 


• 

I 

I 

I 

I 

I 

I 

I 

I I 
I 
I 
I 
I 

I 


o 

o 

•c 


o 

o 


I 

I 

I 

I 

I 

» 

I 

I 

I 

I I 
I < 
I 

I I 
I I 
^ t i 
I ( 

I 


o 

o 

o 


rg 

o o o 

I I I 

u; tii uj 

0 o o 

O O O 

O lA O 

• • • 

^ ^ fg 

1 I I 



(suB.ipey) yoyH3 
-105- 


RANGE (Meters) 
RUN 13 



-O. 


o o o o o o 


o 


o o o o o o 


O O O 0-0 o 


o o o o o o 


o o o o o o 


t/i ^ ^ t 


X 

i 9 

3 

o 


Ui UJ 

o o 
o o 
o o 
o o 


o o 
o o 
o o 
o o 

irv irx 


oooooooooo 


• o o 

tu UJ 

o o 
o o 
o o 
O o 


Of u/ 
o o 
o o 
o o 
o o 


oooooooooo 


X 

o 

o 

u 

o 

a 

o 


o o o o o o 


o o o o o o 


o o o o 
o o o o 
o o o o 
o o o o 
<D « <v N 

oooooooooo 


^sf 

Ui (w U; 

O <N o o 
o o o o 
UJ o o — 

►- o o o o 

< «o m 

X •••••«•••• 

•- oooooooooo 


o 

o 

u 


o 

a 


o 

o 

o 

II 

m 

z 

3 


o o o o o o 

^'<M <V A/ O 

o o o o o o 

Uj UJ UJ UJ Ui UJ 

o o o o <> o 
o m o lo -r lA 

if> ^ tn m 

• ••••• 

o o o o o o 


if) 

*- ^ 
z o 
UI 

T U.' 

o o 
u; c 
«o m 


• • X o 

O'O H- 


a 

*- 

X 


^ IN CM 

o O O O o 

UJ u.' Uj u.' UI 

O O C* lA O 

o o o r* »r» 

IN 00 M ^ ^ 

• • • • • 

o o o o o 


o 

O -.r fv m 
u o o o o 

O K UJ U> Ui UJ 

z o o o o 

•- o o o o 

o o o o o 

Z o O O m 

UJ r- r- -o >o 

oooooooooo 


I UI UJ UJ UJ 
I KJ o Al o O 
O O O o 
O O o o 
o o o o 

' ir» ^ 

I •«•••••••• 

oooooooooo 


u o 

< Z 
U ■“ 


< 

• 

• 

u. 




X 

I I ^ 

o in 


Z 

o 

o 


*- 

K 


m 

in »n 

z 


z 



►- 

or 

or a 


^ u. 



UJ 



u. 

< 

<. 

< 

z 

z z o 

U; 


►- 



<>- 



M 


O O q: 

z 

or 

or 

or 

in 


VI 

u. u. u. or 

UJ < 

UJ 

< UJ 

UJ 

u 





UJ 

u u. 

S> 


a 

s 

O' »- 

►- 

►- 


H- ►- Oj 

u or 

»• 


a 

o 


< 

< 

< 

< 

< < X 

UJ D 3 


3 


< 





n 

rr in 

z 





X 

> KJ 

X 

>> fSl z 




UJ u: u> UJ 
o o o o 
X o o o o 
o o o o 
o o o o 

^ ^ (N O 


— ' oooooooooo 


» ^ ir ^ fw ac o 


- 106 - 


RUN 14 




PLOT OF ANGULAR ERROR VERSUS RANGE 



P.O.OOE 02 6.000E 02 l.OOOF 03 l.^OOE 03 l.BOOE 03 

A.OOOE 02 8.000E 02 1.200E 03 1.600E 03 2.000E 03 



SNUMB = DEDOO. ACTIVITY « = 01. REPORT CODE ='12. RECORD COUNT = 000270 


o 


o o o o o o 


o o o o o o 


o o o o o o 


o o o o o o 


0 o 

1 1 
iLi ui 

o o 
o o 
o o 
o o 

lA 


I I 

4U ui 

o o 
o o 
o o 
o o 


oooooooooo 


rj fsi 

o o o o 


o o o o 
o o o o 
o o o o 
o o o o 
^ ^ ^ — 


^ •••••••••• 

^ oooooooooo 


o o o o o o 


o o o o o o 


o; 

iu 

a 


o o o o 
o o o o 
o o o o 
o o o o 

O) dO Al 

o o o o o o o boo 


a % 
d o 


o o o o o o 


o o o o o o 


Ui (U UJ lu 
O O lA O 

o lA o 

^ lA ^ «M 


0^0 
lA =r «A 
»A fA 


o o o o o o 


^ «A Al =>1’ 4 
Z O O O O ( 
Ui 

z ui ui u; ui 
000^0 
uj o <1 r> o 

«A lA Al a> ^ 


Ui U' 

O. lA O 
O IT* 


XO 0.0 o o o 


Ui 

or 

\J 

z 

XX — 
(O (A 
►H It 

z z o 


u. z 

K H ►- 
»-aaa — 
u. < < < z 


0:(AiOAAU.U.U.Or. 


UJ Ui Ui u 

O M o O 
o o o o 
o o ^ — 
o o o o 

00 lA — — 


— oooooooooo 

O I 

a 

o 

O 04 

o o o o o 

O X UJ U< Uj Ui 

z o o o o 

— o o o o 

O O O o o 

Z O O O «A 

Ui r~ 'O o 

oooooooooo 


^ CA fA ^ 

o o o o 

) Ui Ui Ui Ui 
J IM O <M O O 
O O O O 

C O O O o 
o o o o 

lA ^ — — 

OOOOOOOOOO 


c 

a 

o 

o 

«A \J 

Ui 

cr O 


O tA 

z 


Ui Uj Ui Ui 
O O O O 
X O o o o 
o o o o 
o o o o 

— — AJ ^ 

oocooooooo 


—<<<<« 

< 

X X rsj X 3 


-JUffl— AJfA.tlAOr>40(7‘0 

u. a z ^ 

UiS 3 
a *A z 


- 109 - 


RUN 15 



PLOT OF ANGULAR ERROR VERSUS RANGE 
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plot of angular error versus range 
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COMMENTS ON ILM APPLICATION 


The MADGE equipment as it now exists is not capable of Category III operation. 
Severe multipath environments can cause errors to exceed the limits for 
Category III equipment, focusing errors cause excessive geometric biases, and 
the system is not configured to provide roll out guidance. The following 
changes should be made to MADGE to provide a workable Catetory III equipment: 

o The interferometer arrays should be focused. 

o The azimuth array should be replaced with a space diversity array. 

o The go-around azimuth unit should be mounted back to back with the 
regular azimuth unit, and the pair positioned at the approach end 
of the runway as shown in Figure 22. 

o The arrays should be extended to provide better multipath averaging. 

By mounting the go-around unit back to back with the azimuth unit, its 
measurements can be used for roll out guidance. 

A space diversity interferometer system employing these features can function 

successfully as a Category III_ independent landing monitor in extreme multipath 
environments . 

The coverage of MADGE with these modifications would be unchanged from the 
basic MADGE system except that the extended array would allow elevation angles 
down to zero to be covered, i.e., the coverage is ± 45° azimuth, 0° to 26° in 
elevation, and 50 Km in range. System saturation and mutual interference 
are approximately the same as any airborne controlled DME. The only 
significant difference between MADGE transmissions and normal (VORTAC) DME, 
is the extended reply code containing angle data. MADGE could thus be 
expected to saturate when about 180 aircraft are interrogating one system. 

The equipment as configured has 12 channels to prevent nearby installations 
from interfering. This would be sufficient until large numbers of MADGE units 
became operational. More channels could be allocated if necessary. 
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Figure 22. Location of Equipment for MADGE 
ILM Application 




Some additional processing is required to foc us the interferometer arrays . 
Thus the transponder could be expected to be somewhat larger than the current 
MADGE system. Both elevation and azimuth antennas would have to be larger 
than the basic MADGE to allow Category III operations. The approximate 
parameters of a modified system are given in Table 10. 

The production cost for the complex interferometer /DME concept (Figure 2) 
was estimated by breaking it into it's major components and estimating costs 
of these. Cost of the airborne unit is as follows: 


DME and Antenna 

Data Link Decoding and 
Signal Averaging Logic 

Cross Pointer Display 

TOTAL 


- $ 8 , 000.00 
- $ 4,000.00 


- $ 1 , 000.00 
$13,000.00 


The cost of the ground unit for one azimuth and one elevation antenna is 
as follows: 


20 Receivers at $1,000 ea. 
20 Antennas at $500 ea. 
Data Processing Logic 
Ground Transmitter 
TOTAL COST 


- $20,000.00 
- $10,000.00 

- $ 5,000.00 

- $ 6 , 000.00 
$41,000.00 


A second set of azimuth/elevation antennas, if required for touchdown and 
rollout Category III operations would cost another $35,000.00 or a total of 
$76,000.00 for the ground system. Estimated weight of the airborne unit 
is 8 kgms. 
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Table 10 


MODIFIED MA.DGE PARAMETERS 


Band 

C-Band 

Power 

200 W 

Bandwidth 

3 MHz . 

Size, Weight 


Aircraft Unit 

2-3/8 ATR short 
8 Kg 

Ground Unit 


Azimuth Array 

2 M X .3 M X .8 M 
80 Kg 

Elevation Array 

3 M X .4 M X .3 M 
50 Kg 

Transponder 

,6Mx ,5Mx ,5M 
30 Kg 

Prime Power 


Aircraft Unit 

100 W 

Ground Unit 

350 W 
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ANALYSIS OF A CIRCULAR -SYN THETIC APERT URE LOCATION^SYSTEM-FOR-JLM 


INTRODUCTION 

This study analyzes a circular, synthetic aperture location system called 
"Astrolabe" being developed by the Office National d 'Etudes et de Recherches 

I 

Aerospatiales (ONERA) , Chatillon, France. Although results are derived for 
the capabilities of this specific application, the analysis of the concept 
shows that the results are representative of the much broader class of 
synthetic aperture systems deriving resolution improvement by rotating 
antenna element pairs. The precision of the system is directly proportional 
to antenna element separation and has the potential for providing independent 
landing monitor (ILM) "back-up" or redundancy for low visibility landings 
with the adverse weather conditions and accuracies as specified in Reference 1. 

CONCEPTUAL DESCRIPTION 

The concept by which Astrolabe works is based on the grating lobe structure 
of a two element array where the elements are separated by a distance, R, 
much larger than the wavelength, X, of consideration. A representative 
grating lobe pattern is shown in Figure 23. The pattern is shown in two 
dimensions although it should be considered that a similar pattern exists 
in the third dimension also. The phenomena that allows the Astrolabe system 
to resolve the grating lobe ambiguities is the non-uniform separation (in 
angle) of the grating lobe structure. More specifically, it is shown in 
Appendix C and Reference 5 that the multi-dimensional correlation function 
of the received grating lobe structure can be represented by a Bessel 
function. Appendix C shows the particular matched filter realization of the 
electro-optical system used in the Astrolabe receiving system and demonstrates 
how the apparatus can be altered to modify the correlation function shape 
and trade off correlation function anomalies. 

The correlation processing may be visualized as the movement of a "mask" 
which is representative of the grating-lobe pattern signal being received. 

The movement of the mask to a position which optimizes the registration 
determines the angular coordinates of azimuth and elevation angle. The 
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Figure 23. Antenna Grating Lobe Pattern 
Representation 
. -117- 



_analogy of a_^mask!LJ.s_used_because— the eleGt-ro-opfeieal— apparatus~in— the 

Astrolabe system uses an optical transparency mask to perform the correlation 
by synchronizing the movement of the mask to the spinning of the antenna, 
uniformly illuminating the mask with a light source with intensity proportional 
to the square-law detected signal, and observing the correlation pattern 
produced as the eye Integrates the light transmitted through the spinning 
mask. 

The concept is valid for either a spinning transmitting or receiving antenna 
pair. This can be seen from the expressions in Appendix C for the two signals 
(one from each antenna) input to the detector. 

A^ exp j t - ^ ) 

c 

A 2 exp j (fi t - ^ ) 

c 

The signals depend only on the absolute separations, d^^ and between the 
transmitting and receiving elements. Similarly, doppler shift expressions 
(which modify f2) depend only on relative motion between elements. The 
system concept is thus unchanged whether a single, fixed-gain receiving 
antenna detects a revolving, grating-lobe transmitting pattern or a single, 
fixed gain transmitting antenna is detected by a revolving, grating -lobe 
receiving antenna pattern. 

GEOMETRICAL CONSIDERATIONS 

The basic system concept described in Appendix C obtains angular Information 
by interpreting the received signal produced by square-law detection of the 
superimposed signals from two revolving antenna elements. The geometries to 
be considered may consist of rotating transmitter or receiver antenna element 
pairs either on the ground or in the aircraft. The angular Information is 
available only at the receiver unless a separate data link is provided. 

The angular information consists of azimuth, 0^, measured in the plane of 
rotation of the antenna pair (with an arbitrarily chosen "zero" reference) 
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and elevation angle, relative to the plane of rotation (plane of rotation 
is defined as the plane normal to the axis of rotation and containing the 
center of the array.) 


Background 


Reference 6 derives the expressions for the correlation function and relates 

each term of the expression to the correlation phenomena produced in the 

electro-optical receiving and display apparatus. The azimuth and elevation 

angle determination is made on the resolution of points with an intensity 

profile described by the Bessel function and + g.2]^ 

The two terms correspond to a false target at - 0 g and and the true location 

at 0 , ip . Assuming that the system is initialized to remove the false 
s s 

correlation point as per reference 6 (i.e., a + g = ir/2), we have as the 
expression for the true location point intensity profile: (See Appendix C 

for definition of terms . ) 


C J + b« 2 ] 

0 0 


where + 3*2 _ 


ir^p2 


Stt^PR , ^ . 16tt^R^ 2 . 

cos Ip cos (0 -0 ) + = — cos Ip 

6 \ s s o ,2 s 


2 „ X 16ir^R^ 2 . ^ 2 „ 

cos (0 -0 ) H 7 cos ip sin (0 -0 

so X ^s so 


= 1 ' 

_ 


Making the substitution P = —7 — cos i|;. 


IT p 8 tt PR , /n r. \ . 16 tt R 2 , I h 

^ ;t— cos |() cos (0 -0 ) H r COS \p ( 2 ) 

-2 S so ,2 s' 


46 R 




+ B'^ = [cos^ ip - 2 cos ip cos (0 -0 ) cos \p + cos^ ip ]'^ 


s o 


Ip and 0 are the desired angular references (actual values) . 0 ^=wt + 0 

s s ° ant s 

is the time dependent reference of the spinning antenna. The disk which 
performs the cross-correlation is spinning with the same angular frequency, 
w, such that its angular time expression is + 0 ^. The relationship 

between 0^ and 0^ is thus some constant which allows expression of angular 
position, 0 ^, on the disk to the true azimuthal position 0 ^. The value 
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/ 


6^-6^ in the above expression is a cons tant value whlch_ls-chosen-^to-slmpl-tfy 
calculations. 


Reference 6 provides a description of the coefficients of the Bessel function. 

r~i — T 

/A' + B' are described as vector magnitudes. 


The magnitude 




These coefficients, + B'^ and 


2 

+ B' , is representated as |u - v| where v describes the 


receiver location and u describes the observed correlation point on the disk. 

These vectors are shown in Figure 24. To analyze the effects of spot size 

on the disk as a function of angular position, the vectors are set equal, 

u = V, such that 0 =0 and |u - v| approaches 0 for tj; approaching ij; . The 

s o y 2 ® 

previously derived expression for </k' + B' 


becomes (for 0 -0 =0) : 
s o 


y/k'^ + B'^ = 


4irR 


(cos \b - cos ) 
s 


(4) 


are e qual, u = v, 
J i/k'^ + B'^] a 


the spot profile is independent of 0 and 


as ijj approaches ijj is descriptiv e of the 

+ B'^ = 1.55 


When the vectors are e< jual, 
the variation of 

spot intensity variation with elevation angle. Since 
corresponds to the 3 dB point of the Bessel function, this definition of 
resolution may be used to derive the value of \p corresponding to the 3 dB 
point and thus derive the precision of the correlation function. 


Thus, 


1.55 


= + B'^ = 


or <1/ 


3dB 


cos 




4ttR 

X 

55 X 
4ttR 


(cos ij) - cos \p ) 
Job s 


+ cos 


5 


(5) 


( 6 ) 


Subtracting ij; we obtain the angular uncertainty corresponding to the 3 dB 
s 

point of the correlation function. 


AiP = - iPg (7 

The relationship of Aip as a function of is shown in Figure 25. 

The precision of the correlation function in the 6, or azimuth position may 
be determined from the geometrical symmetry of the spot at distance P from 
the center of the disk. From the derivation of the correlation function in 
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Appendix C, it is known that the spot at the point (P,6^) is a circle. The 
dimension in the A\() direction is the same as the dimension in the A0 direction. 
The radius at the spot is : 


distance = 


46R 

A 

46 R 
A 

46R 


cos (ip + Alp) 

(cos iji - A sin ip) 
46R 


cos \j> - 


A A 

The circumference at this elevation is 


Alp sin ip 


( 8 ) 


Thus 


2ir* Radius 


8tt6R 

A 


cos Ip 


2tt 


46R 

= _J. 

8tt6R 

A 


Alp sin ip 


cos ip 


Alp sin ip 
2ir cos ip 


(9) 


( 10 ) 


A0 cos i|; = Aij; sin ip (11) 

This relationship for the azimuthal uncertainty as a function of elevation 
angle and elevation uncertainty is also plotted in Figure 25 to show how 
the uncertainties relate. 


Accuracies and Precision 

Angular precision has been derived in the previous section to describe the 
repeatability of the correlation function spot position. The application of 
the electro-optical correlation apparatus described for the Astrolabe 
reference provides correlation function analysis and interpretation by 
human observer techniques. Operational tests of the Astrolabe concept 
using this apparatus for display and parameter retrieval have demonstrated 
2a accuracies comparable to the anticipated theoretical precision (Reference 5) 
The theoretical precisions derived in the previous section will thus be used 
as a baseline to analyze conceptual capabilities. 


The study of the correlation function parameters identified angular deviations, 
A0 and Aip, which are descriptive of the conceptual ability for measuring 
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'direcITionrwithrrespect^to^the^spinning^antenna^refereiice-frame^ — The-receiver — 
processing is limited to deriving angular information referenced to the 
spinning antenna at the receiver. If the spinning antenna is the receiving 
antenna, the angles may be directly referenced to the receiver's reference 
frame. Otherwise, the angles derived are relative to the spinning transmitting 
antenna and additional reference angles may be required in order to orient the 
two separate reference frames. Table 11 provides a comparison of the different 
geometrical configurations for concept implementation. Systems consisting of 
one or more antenna placement geometries will derive guidance or position 
fixing by three techniques; location corridor, corridor triangulation, and 
geometrical shape orientation. 

Location Corridor 


The Location Corridor for monitoring consists of establishing directions, 

6 and tj;, within the angular uncertainties of A6 and Ai);, Linear uncertainties 
are SR*A6 and SR'Atp where SR is the unknown slant range between the receiving 
and transmitting antennas. 0 and il; are referenced to whichever antenna is 
spinning. Figure 26 shows the linear uncertainty as a function of antenna 
element separation and range. 

Corridor Triangulation 

Corridor Triangulation establishes location by intersection of two or more 
independent (separated) location corridors. The corridors defined by 
SR£A6^ and SR^Ai|;^ and SR2*A02 and SR2*Ai(;2 may or may not intersect for each 
set of direction angles measured, AS^^, Aip^^, and A62, Ai|^2* Assuming that 
sufficient integration time is available to provide corridor Intersection 
as defined by the precision A 0 , A\|; and eliminating non-intersections as 
false multipath targets, the volume defined by the Intersection will have 
a maximum cross-sectional linear dimension of the larger of: 


and 


SR-Aij; 


SR-A 9 

, . -1 D . 
cos (sin 2^) 

2 -SR^-A 6 

D 


A angles in radians 
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Table 11. Geometrical Configuration Consideration 
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LINEAR UNCERTAINTY, METERS 
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Figure 26. Linear Uncertainties for A 0 = A\l> (i|j= 45°) 
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These expressions are derived from Figure 27. An independent landing monitor 

will require resolving angular deviations in both horizontal position and 

altitude if dependence on other systems (such as a radar altimeter) are to 

be avoided. Equal azimuthal and elevation angle precisions will be assumed 

to provide comparable accuracies. Altitude accuracies (two a values) within 

the SR- Alp precision should be expected with similar values for azimuth 

accuracy (See Figure 26) at longer ranges (2*SR>>D) degrading as dj^. Range 

accuracy for this assumed configuration with the location corridor sites 

straddling the runway would thus vary as d_. Figure 28 is a graph of the 
d ^ 

factor which relates range precision to azimuth precision (and similarly, 
*^1 

the accuracies) so that expected ranging performance may be anticipated by 
using the graph of Figure 26 multiplied by the factor in Figure 28. 


Geometrical Shape Orientation 

The third possible technique for guidance or position fixing consists of 
determining aircraft location by measuring relative directions to a known 
geometry. An example of this woixld be beacons defining a runway edge. The 
two primary geometrical dimensions used for defining azimuthal (and elevation) 
angle deviation and range are shown in Figure 29 . The angular deviations 
derived as a function of range will assxane that point pairs used for range/ 
roll attitude and azimuth/elevation angle derivation are known so that 90° 
approach ambiguities are not possible. Although the geometry is idealized, 
the results are indicative of guidance and positional capabilities. 


The dimensions L and W represent the geometrical point separations . When 
0 >>A0 (which may be assured by selecting W sufficiently large for the 

required detection ranges). 



then Range = 


W ± 2«Range»A9 
2 • Range 


W 

y + Range* A0 



( 12 ) 

(13) 
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Fig. 27. Location Corridor Site Geometry for Corridor Intersection 
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Figure 28, Relationship of Range and 
Azimuth Precision 
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GEOMETRY PLAN VIEW 



Figure 29. Geometry for Deviation of Guidance and 
Position from Known Ground Points 



or Range 


Range 


1 + A6 


sin ^ ^rangej 


- W 


e 


2 sin range 
2 


1 + 


2 -Range A0 
W 


2 


sin I range \ 

i ^ ) 


(14) 


(15) 


This expression may be evaluated to find the two ranges satisfying the limits 
The total linear range variation, A Range, as a function of the actual range 
and the separation, W, may be defined as the difference of these two range 
values or: 


» D 4 -Range A0 

A Range 2 ^ 


The angular deviation from runway centerline, AAz is; 

-1 


AAZ 3 sin 


L cos ijj 

(Range + 2 ) A0 


L cos ij) 

The elevation angle deviation, A El is derived below. 


and 


I . = (1 

min 


L 2 2 

(Range -cos - y) 4- (Range -sin \p) 


R 


max 


L 2 2 

= (Range. cos + y) + (Range -sin ip) 


U. 


0 ' 

■]' 


(16) 


(17) 


(18) 


(19) 


For each of these ranges define an elevation angle ij; as shown on Figure 29b. 
,-l Range -cos if; -L/2 


= cos 


ip2 = cos 


min 

-1 Range -cos i); +L/2 


R 


( 20 ) 


( 21 ) 


max 


The deviation of the point locations within the envelopes of Atp may vary the 
derived value of elevation angle, about the true value, ip , such that: 


♦der 2 * 


_1 R I 

^ tan I min- Alp _j_ max - Aip I 

1_L cos ip^ L cos ip 2 J 


( 22 ) 
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= ± 


tan 





fR . 

R . . \ 


Ml 

mxn 

1 max 1 



i cos 

cos Jj 

(23) 


The total geometrical angular deviation is thus 


Ae = tan 


-1 


min* Aiji 
L cos ij) 


max • Alp 
L cos 1 ^) 2 ] 


(24) 


Figure 30 shows the factor that is used to extend the results of Figure 26. 

It should be noted that the factor A Range is approximately twice the factor 
in Figure 28 such that range precision and accuracy of ranging for geometrical 
shape orientation is worst (approximately twice as large as location corridor 
intersection) for the same goemetrical separation (W versus D) . 

Figure 3l shows AAz, the angular deviation from runway center line for 3° and 
6° glide slope approaches. Guidance point separations, L, of 8 Km, 4 Km, 
and 2 Km are shown. The values are independent of ip for i|) < 6°. The error, 
or range of deviations, are ±AAz. For small angles, AAZ<10°, AAz is directly 
proportional to A0, the azimuthal angle precision of the correlation function. 


Figure 32 shows the elevation deviation. The curves are approximated in 
the region of sharp upsweep caused by overflight of the nearest geometrical 
reference point. The curves are insensitive to variations in \p for t|;<10° 
and are directly proportional to Ai|; for Ai|k 1.0°. It becomes apparent that 
better elevation angle precision (almost by a factor of two) and a series of 
geometrical reference points are required to achieve the precision of a 
corridor location technique for glide slope guidance. 


Conceptual Implementations 

The applications for guidance or positioning by the three methods of corridor 
location, location corridor triangulation, and geometrical shape orientation 
are reviewed in Table 12. One conceptual implementation to achieve A0 and 
AtJ; angular precisions relative to each spinning antenna has been described • 
as the Astrolabe system. Other implementations of this circular, synthetic 
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14 . 



Figure 32, Elevation Angular Deviation as a 

Function of Geometry Point Separation 



aperture concept would not differ substantially as far as transmitter, 
receiver, and antennas are concerned. The key parameter of the system would 
still be an antenna element separation of several wavelengths in a multi- 
element antenna which generates a multi-lobed antenna pattern which is 
rotated, or moved, in a predictable manner. 

The number of lobes resulting from the multiple wavelength separation is 
directly a function of the separation distance so, similarly, the correlation 
function periodicity (and indirectly, the 3 dB point) are more a function 
of physical separation than individual element type or pattern. The Bessel 
function correlation characteristic is derived for a specific receiver 
processing application, however, it is apparent from Appendix C that other 
matched filter implementations would provide similar correlation characteristics 
(precision) differing only in relative strengths of the rings surrounding 
the correlation point, suppression of the ambiguous point, and biasing as 
a function of elevation angle (the same factors affecting the electro-optical 
apparatus). Spot size and ring diameter, which are related to precision, 
would be essentially independent of the processing scheme used. 

The electro-optical processing described in the Astrolabe documentation is 
the most economical of several processing alternatives. The major objection 
to this application is the requirement for a trained observer to interpret 
the visual correlation function bo derive necessary direction angles (location 
corridor descriptors of azimuth and elevation angle) . This task involves 
rejection of multipath target angles which are unresolvable without knowledge 
of their time history. The unpredictability of this multi-path mechanisms 
makes automated processors difficult to implement due to the complex pattern 
recognition task of point pattern interpretation. Reference 5 describes 
how cross-terms (off-diagonal) of the multipath signal matrix are resolved 
by time integration, however the direct multipath terms are unresolvable in 
the same integration interval. 

These multipath terms must be resolved from the true direction angles if 
precision guidance or positioning is to be achieved. Within the scope of this 
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-analysis— there— is insufflcient-knowledge of—the~characteristics^of“the“"direct 
multipath interference to specify the requirements of an automated decision 
processor. It can be assumed that such processing would not be more effective 
than a trained observer since the processing would be against the same 
correlation function as that viewed by the observer. For these reasons, the 
human observer/interpreter is proposed for the ILM application with an 
automated system as an alternative with greatly increased cost and complexity 
with no increase in system precision. 

ILM APPLICATIONS 

The concepts used in the Astrolabe system have potential for use as an ILM. 
Reference 1 contains the desired ILM angular accuracy specification of 0.10° 
to 0.12° for 2a accuracy. Reference 5 presents some experimental results 
of measurement accuracy of a prototype Astrolabe system wiLh an R/X=20 
normalized element separation. The empirical values for 2a accuracy were 
within approximately 0.4° for azimuth and elevation angle determination. 

Relating these to the theoretical precision of 0.50° for R/X=20, we can 
expect the 2a accuracies of the application to be within the predicted precision 
for the respective R/X value. To assure 2a accuracy of 0.12° would require 
element separations of approximately R/X=80 which would require spinning 
element separations of approximately 3 meters at X-band and 0.8 meters at 
K^-band. The selection of frequency for the transmissions is somewhat 
arbitrary although it should be recalled from Reference 1 that approximately 
(1.3 X R (Km)) dB more power is required for K -band operation than X-band 

IHclX d 

operation in heavy rain (16 mm/hr) to compensate for path attenuation effects. 

For example, guidance or position determination at a range of 8 Km would 

require over 10 dB (10 times) more power at K -band than at X-band for the 

cL 

same receiver signal to thermal noise ratio. 

It is conceivable, however, since an 80 cm element separation would be possible 
in the nose of a commercial aircraft, that K -band operation of a spinning 
antenna in the aircraft for high accuracy angle determination would be feasible 
for configuration 1.B.2 or 2. A. 2 of Table 12. Both these configurations 
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would have the requirement for stabilization of the spinning antenna to 
eliminate high frequency roll, pitch, and yaw angular deviations to achieve 
the required high accuracies. A gimballed system for the spinning antenna 
(with the spinning elements generating the gyroscopic stability platform) 
could be used with sensors to derive smooth aircraft to antenna reference 
angle to relate derived angles to the aircraft reference frame. 

Several references have been made previously to multipath errors. These 
consist of two types: direct multipath terms which are unresolvable from 

true targets over short observation times and cross-product multipath terms 
that originate from the product detector which are resolvable over relatively 
short (less than one second at X-band) integration times. Reference 5 contains 
a detailed analysis of the effects on the correlation function for high 
reflection coefficients. In summary, the cross-product terms are reducible 
by at least 12 dB and the direct multipath terms are left to the expertise 
of the display analyst for removal. The multipath problem therefore favors 
ground receivers and analysts since the multipath geometry would be more 
familiar and predictable. 

Comments on ILM applications up to this point have been general for nearly 
any configuration. When multiple systems are used in an effort to attain 
location by triangulation or geometrical shape orientation, several other 
considerations are necessary. 

Multiple ground receiving locations have been shown to require large separa- 
tions (2 to 4 Km) . Range factors (the ratio of range uncertainty to linear 
azimuth uncertainty) are to be kept below 10 for ranges up to 10 Km. 

Reference 1 specifications identify required range factors of 2 to 3. These 
wide separations seriously hinder displaced ground receiving sites (configura- 
tions 1.A.2 and 1.B.2 of Table 12) since the triangulation function must be 
performed centrally. Either multiple receivers, receiver processors, and 
analysts or multiple receivers with IF or detected IF signals relayed to a 
central processing and analysis center would be required which would be much 
more costly than some of the alternatives to be presented below. 
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transmitters have the same problem of wide separation of ground antennas for 
high accuracy although the analysis is centralized to the aircraft. The wide 
separations necessary for ranging accuracies necessarily limit the linear 
accuracies to D/2 • A0 , D/2 • Aiji, W/2 • A8, and W/2 • Aiji. This effect coupled 
with multipath problems will degrade performance such than any multiple 
configuration will be of questionable use for flare, touch-down, and roll- 
out information. 

The multiple ground antenna configurations, although capable of providing 
high accuracy results prior to flare maneuvers is not an economical ILM 
configuration because of the complexity of the analysis, interpretation, and 
triangulation requirements. If the high accuracy requirements can be 
waived it would be feasible to use configuration 2. A. 2 where the imagery (4 
beacon points) would be representative of the four corners of the runway and 
presented to the aircrew in the form of a real-world perspective image for 
a VFR type of approach. This application would require the spinning (or a 
stepped Wollenweber array) receiving antenna in the aircraft. Decreased 
accuracy requirements to the order of 1 degree would negate the requirement for 
the stabilized antenna platform. Operation at lower frequencies and/or 
smaller element separations would also be possible. 

A high accuracy ILM application at a reduced cost when compared to the multiple 
ground antenna configurations would consist of a combination of a single 
beacon transmitter on the aircraft, a spinning receiving antenna on the ground, 
and an interrogator /response DME (Distance Measuring Equipment) ranging 
capability incorporated into the ground-to-aircraf t data link for relaying 
angular information to the aircraft (Configuration l.A.l extension). This 
system would provide high precision approach corridor information. (Azimuth 
and elevation angle along with precise range- to-go as described in Reference 1, 
ILM Multilateration Sensor section.) It should be pointed out that this 
combination of sensors would benefit from the most accurate capabilities of 
both techniques: Astrolabe precision angle referencing without triangulation 

techniques and its associated accuracy degradation and high precision ranging 
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by TOA (Time-of-Arrival) techniques without GDOP (Geometrical Dilution of 
Precision) inherent in multilateration angle determination. 

Configuration 2.A.2 is a less precise system than Configuration l.A.l 
extension and provides a pseudo real-world perspective image for a visual 
back-up ILM. Point accuracies of the order of 1° allow a qualitative four- 
point image for flight crew reference. Configuration l.A.l extension is a 
high precision (2a angular accuracies of approximately 0.1°) angular Astrolabe 
array with DME ranging accuracies of the order of a few meters. 

The physical equipment requirements and functions of the two variations are 
summarized in Table 13. 


Mutual interference of variation 1 is non existent when the beacons are 
commutated "on-off" in a cyclical pattern. Secondly, all the airborne 
equipment is passive (receiving only) so the same runway beacons may be used 
by all aircraft without conflict. Variation 2 would use coded airborne 
beacons with the ground station "repeating" or re-transmitting the code 
identifier so the appropriate aircraft may receive its coordinates. The only 
mutual interference would be in the multiple usage of the data link/DME data 
bus which would have to be maintained and controlled by the ground station to 
assure user coverage. This would be possible by using angle and range gating 
techniques. For example, a DME measurement exchange could be used to turn 
beacons on or off as required to derive angular coordinates and give update 
rate priorities to aircraft nearing touchdown. 


Projected costs of these systems were generated by estimating major component 
costs. For variation one (Figure 3), the airborne costs are as follows: 


Receiver 

Antenna 

Signal Processing 
and Display 

TOTAL 


$ 2 , 000.00 
$ 5,000.00 
$ 4,000.00 

$11,000.00 
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Table 13. Equipment Svnnmary for ILM Applications 


Variation 1 Cockpit Perspective Display 
Ground Equipment 

4 beacon transmitters with antennas 
Commutated transmission control 

Airborne Equipment 

Rotating array 
Receiver 

Electro-optical processor 
Analyst 


(Configuration 2. A. 2) 

Usage 

Outline runway 

Coordinate beacon transmissions 
Usage 

Provide synthetic aperture 
Detect beacon signals 
Display correlation function 
Interpret position 

(Configuration l.A.l extension) 


Variation 2 Ground Station Analysis and Control 


Ground Equipment 

Rotating array 
Receiver 

Electro-Optical processor 
Analyst 

Data link/DME 

Airborne Equipment 
Data link/DME 

Ins trumentation 

Beacon transmitter with antenna 


Usage 

Provide synthetic aperture 
Detect beacon signal 
Display correlation function 
Interpret position and enter to 
data link/DME 

Measures range and transmits 
position information to air- 
craft 

. Usage 

Measures range and receives 
position coordinates from ground 
Displays aircraft coordinates, 
range-to-go, and/or Az-El 
angles from runway 
Provides signal for ground 
fixation 
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This concept Includes four ground beacons at each corner of the runway. The 
cost of each beacon was estimated at $5,000 for a total ground system cost 
of $20,000. It should be noted signal processing for this system consists of 
a mechanically spinning reference transparency with a light source driven by 
the received RF signal and a simple display of the cross-correlation function. 
This probably is not adequate for an operational configuration. 

Variation two (Figure 4) consisted of an airborne beacon transmitter with 
ground receiver, processor/display and a DME/data link for ranging and transfer 
of guidance data back to the aircraft. Costs for the airborne unit for this 
variation were as follows; 


Beacon - $ 6,000.00 

DME/Data Link - $ 8,000.00 

Cross Pointer Display - $ 1,000.00 

TOTAL $15,000.00 


Ground system costs are as follows: 


Receiver 

Signal Processing/ 
Display 

Rotating Antenna 

DME Transponder and 
Data Link Logic 

TOTAL 


- $ 4 , 000.00 

- $ 4,000.00 

- $ 5,000.00 

- $10,000.00 

$23,000.00 


Estimated weight of the airborne unit for this variation is 8 kgm. This 
variation also has the low cost mechanical signal processor at the ground 
site which probably would not be suitable for an operational configuration. 
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Section V 
Appendices 
Appendix A 

Determination of Array Gain Characteristics 

The generalized expression for the distant field of each element of an array is 
(See Figure A-1) 

Ei(0.<t')=f(e,'|))li exp p(kri cosii^t + 0 ^ 4 ^ (A-i) 

with: 

- the far-field function associated with the individual antenna element 
It ~ amplitude excitation 
~ phase excitation 
k s 2ir/X 

cosipc = COS0 cos0i + sin0 sin0i cos(<j>-<j)i ) 

r^ = range 

The total field contributed from an array is the summation of the individual 
contributions of the n elements. 


n 

E(0,<j)) = Z Ei(9,<l') 
i = l ■ 

= f(0,<fi) z It exp lj(kri cosYt 
1 = 1 ^ 

= f(e,<p) • S (a- 2 ) 

f( 0 ,(fi) - element factor 
S - array factor 


The total field may be considered as consisting of the two separable factors, 
the element factor, f(0,<f>) and the array factor, S. The power gain of the array 
may be expressed as the ratio of the squared field intensity in a particular 


direction to the average radiated power intensity. 

G(e.*) = , 4. | f(9,4.) | ^ | S | ^ ~ 

P/4it P 


Where P is the total power 
delivered to the array 
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The element factor of a typical array element applicable to an ILM 
aircraft mounting is presented in Reference 2. A linear approximation 
of the vertical and horizontal gain curves presented in Reference 2 
are shown in Figure A-2 and A-3. These curves were numerically 
integrated to derive the lossless antenna maximum element gain as 
20 dB. The element factor, |f(^)[ has the normalized gain character- 
istics shown in Figures A-2 and A-3 with a maximum gain of 20 dB at 
boresight (fmax ^max additional 

gain due to the array factor, S, will be analytically derived in the 
following sections. 

The concept of a thinned array dictates that elements are removed or 
"thinned" from a conventional array with element spacings less than 
the wavelength, Removal of these elements creates grating lobes. 

The uniform amplitude distribution, = 1 has the unfortunate property 
of high sidelobes (only 13.2 dB below the main beam). A bistatic thinned 
array removes the grating lobes to null, or at least lower sidelobe, 
regions of the transmitting antenna, however, the receiving array with a 
beamwidth much narrower than the transmitting antenna, must sacrifice some 
mainlobe gain (and beamwidth) to reduce the sidelobes to an acceptable level. 
The amount of reduction is somewhat arbitrary, however, it can be assvimed 
that if good runway edge resolution is required the sidelobes must be below 
expected edge contrast levels. The runway scenario of the previous studies of 
imaging radar ILM's (Reference 1) identified typical concrete runway to grass 
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Figure A-3. Normalized Antenna Gain Profile (Elevation) Elevation Profile 
(Reference to wing surface in the direction of flight, 4-angles 



contrast ratios at X-Band in excess of 20 dB. 


As suming a t— iea s t^6-dB^addit i onal— different iat^fo r^de t ec t ion~c r t ter ia^ 

the sidelobes must be at least 26 dB dovm from the main beam. The 
proposed array in Reference 2 is in close agreement with Taylor weighting 
on the anplitudes for sidelobe suppression to approximately 33 dB 
below the mainlobe. The effect of this amplitude weighting is to produce 
a slight decrease in the maximum theoretical mainlobe gain of 2 to 
3 dB. 


This factor will be incorporated later. 

g ^ max lE(e,(p)l 

\ax P/'4n 

4n max | E(e,(j)) ( 


2n n 2 

/ / 1e( 6.<I>)1 sine de d(j) 


0 0 


(A-4) 


or G 




2 ,2 

4n max |f | |$| 


ax 2n n 22 

/ / |S| sine de 

0 0 

Previously we derived the maximum value of the element gain by: 


(A-5) 


max 


_ max |E(e,(^)i 
P/4n 

4n max If (e, 4 >) 


2 n n 2 

/ / |f(e,<j>) 1 sine de d<}. 

0 0 


20 dB 


(A-6) 


for |Sp = 1 (one element only) 
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This expression was evaluated on the assumption that 






G 


max 


max |f(<}))| 

= 2n 2 

/ |f(<(>)| d4> 

0 


max |f(e)r 

n 2 

I |f{e)l sine de 


= 4 n 


G ( 0 ) 
max'^' 




(A-7) 


When the element is used as an array element, the array factor must 
be derived. This gain will have a maximum value expressed by: 



[ 2 2 

|f(e,({)) I |s| ] 

2 n n 2 2 

/ / lf(e,<t>)j Isj sine de d(|) 

0 0 


(A-8.) 


When f( 0,<ti ) and S are maximum in the same direction: 


Ga, 


max 


4 n max|f(9,(j)) 1^ maxjSl^ 


A-9) 


Previously it was described how the pattern of the array has narrow 
beamwidths with several grating lobes, the number, m, of lobes being 
determined by the element spacing. The array factor effectively shapes 
the array pattern since the array bearawidth in 0 is very narrow. 

Thus, for each lobe the variation of f(0) across the bearawidth is small 
such that (since S 'is independent of 9 for the array geometry chosen) : 


0 de d(j. = /'|f(e)|''ls|‘' sine de d<j> 


I 2 , 7 T I 


Continued 
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(A~10) 



where 


fi(e) 


is the element factor constant for the i 


th 


lobe 


0 i 



is the angle corresponding to the i^h lobe 
(0 < < it) 

is an angle chosen sufficiently larger than the respective 
lobe beamwidth such that neglibible power is radiated beyond 
this angle until the next lobe is encountered. 


Then 




■max 


= 47T 


rs O 

maxlf(4)) r max|f(6)r max|S 






(A 


1 1 


Recall that; 

p 

max 


. 4^ m aK|f(e)|^ 

f lf(e)rsin0 d0 d(j) 


0 


So 


^Amax 


= G 


maxjSl^ lf(6)l sin0 d0 


(i 


max 




The lobes are very similar over the range of the f(0) function when 
f(0) has a narrow beamwidth (less than 40° to 60°) such that: 


^Amax 


= G. 


max 


max S 

L. 

/g®^^^|S|^sin0 d0 

r fi(e) ^ , |f2(e)|^ ^ 

i 

f .. ••• 

- 0 0 


Dividing through by the term f-j (9 ) | _ 6_,^(0) 

max 

^7T ... 

0 


- 11 ) 


.- 12 ) 


.-13) 


(A-14) 
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(A-15) 


’Amax 


max 




max I : 

Qsi' 


sine de 


G-| (e) 62^61 

G (e)”^ 

maw ' ' mav ' * 


max 


max 


■ *"G [F) 

max' ' 


where (0) is the gain associated with the i*"^ lobe. 

When the array is thinned such that the element spacings, although uniform, 
are in excess of one wavelength, the summation becomes the value K' and the 
expression for becomes : 




max 


6 

max 


1 

K* 


max 


sine de 

e-c 


^max K" 


(A-16) 


The values of K' are only a function of element spacing since it is this 
characteristic that generates and places the grating lobes. Values of K' 
as a function of element spacing for the array element azimuth characteristic 
of Figure A-2 are shown in Figure A-4. 

The factor K", with the denominator integral limited to one lobe is the same 

as the integral expression for an unthinned array of length % with no grating 

lobes. Max |s | ^ is the normalizing value associated with the |s| ^ term in 

the integral. K" is thus indistinguishable from: 

o 

max Is I 

* . 2 for a uniformly spaced array with no grating 

|S| sine de 

lobes which has been approximated in Reference 8 as : 
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Figure A-4. K' as a Function of Element Spacing 
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4tt max|S|^ ^ ^ 

2tt /^|S1^ sine de X 
0 ' ' 


for long arrays which is independent of 
the number of elements 


Thus, K" 


= ifor 


Jt » X 


(A-17) 


The expression for maximum array gain becomes: 


^^max 





X 

K'X 


(A-18) 


For the array element shown in Figures A-2 and A-3 and assuming a 3 dB loss 
associated with amplitude weighting as per the previous discussion: 


^Amax dB ' [x] " 

(A-19) 
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APPENDIX B 


A COMPUTER PROGRAM FOR THE 
ANALYSIS OF SPECULAR MULTIPATH IN AN INTERFEROMETER 


I. ENVIRONMENT 

The environment in which the interferometer is assumed to be working 
is a quasi-two dimensional. The antennas are positioned along the Z 
axis , and the aircraft flight path is in the X-Z plane. A flight 
path component is also allowed in the Y dimension. All reflecting 
surfaces are assumed to be finite line segments in the X-Z plane with 
infinite extent in the Y direction, and orientation perpendicular to 
the X-Z plane. The multipath geometry is as shown in Figure B-1. 


The method of images is used to compute the reflection effects. For 
each antenna and each surface, the surface is extended to infinity and 
the image of that antenna about that surface is formed. A line is drawn 
from the aircraft to the image. If the line intersects the extended 
surface within the bounds of the real surface, and the extended surface 
does not separate the real antenna from the aircraft, a reflection is 
presumed to occur. The angle of incidence between the line and the 
surface is then computed, and the reflection coefficient computed from 




COS0 




® g-(2KaCOS0)^ 


COS0 + ’yJe^-sin^Q 


(B-1) 


where 


■^f 


2 

COS0 --\ /e -sin 0 „ 

-(2KOCOS0) 

Rh = e 


COS0 +~\/ep-sin 0 


(B-2) 


is the complex reflection coefficient for vertical polarization. 

Rjj is the complex reflection coefficient for horizontal polarization. 
0 is the incidence angle 

is the complex permittivity of the surface 
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AIRCRAFT POSITION 
(XAC, ZAC, YAC) 



Figure B-1. Multipath Geometry r 

(lypical for one antenna pair, one surface, one aircraft, position)pDai^i&jbositr^ 


a is the RMS 'surface roughness, 

K is the wave number. 

The magnitude of the reflection coefficient is modified by the ratio 
of the range from the aircraft to the real antenna to the range from 

the aircraft to the image, and by an approximation to the antenna 

pattern. The phase of the reflection is modified by the difference 
in these ranges . 

The reflection components for all surfaces are summed vectorially 
with the direct component to obtain a total field present at a given 
antenna. The direct vector is used as a magnitude and phase reference 
as shown in Reference 1, so that the normalized vector sum is a measure 
of the error in phase at the antenna. This is converted to a normalized 
range difference error by dividing by the wave number. 


An estimate is then formed of the range difference between antennas of 
a pair using either a focused or unfocused algorithm. If the run is 
unfocused, the difference between range errors is added to the actual 
range difference. This range difference is divided by the distance 
between antennas to form an estimate of the sine of the angle between 
the aircraft and the x-axis . If the run is focused, the algorithm is 

0 = R^(2)-r2(1) + ZA(1)+ZA(2) 

2R(1) (ZA(1)-ZA(2)) 2R(1) 


(B-3) 


where 

^ is the estimate of the sine of the angle. 

R(i) is the range from the aircraft to the ith antenna. 
ZA(i) is the height of the ith antenna from the x-axis. 
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The estimates for the various pairs of antennas in the array are 


averaged in an unweighted average. Since the sine of the angle rather 
than the angle is available, the error is computed via the approximation 


a; 

e 

Of 


6 - 0-6 - sin sin (©-“S) >"8^0008^-3 inOcosO 


(sin0-sin0) cos© for small 0 
where 0 is the angle between the line from the aircraft to the 
origin and the x-axis . 

0 is the estimate of that angle. 

^ is the error in that angle measurement. 


These errors are computed at each increment along the flight path, and 
output both in tabular form and a plot. 


) 




(B-4) 

(B-5) 
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GL05SAPY OF TEPMS: PROGRAM MAnGE 


RGH 

REAL 

surface RMS ROUGHNESS DIVIDED BY WAVELENTH 

RHO 

REAL 

MAGNITUDE OF FRESNEL REFLECTION COEFFICIENT 

PHOP 

PEAL 

MAGNITUDE OF REFLECTION COEFFICIENT MODIFIED BY RANGE 



AND ANTENNA PATTERN RATIOS 

RR 

DOOPLE 

RANGE FROM AIRCRAFT TO IMAGE ANTENNA 

RRFF 

OOUPLE 

RANGE FROM AIRCRAFT TO REFFRENCF POINT 

SKP 

REAL 


SMEAN 

real 

sample mean 

SMSO 

REAL 

sample mean square VALUE 

SNTH 

OOUPLE 

SINE OF TRUE PHYSICAL ANGLE 

SURF I 

REAL 

INTERCEPT OF SURFACE WITH Z AXIS 

SURFS 

REAL 

SLOPE OF SURFACE 

SVAP 

REAL 

SAMPLE VARIANCE 

TEMPI 

REAL 

TEMPORARY INTERMEDIATE VARIABLE 

TFMP2 

PEAL 

TEMPORARY INTERMEDIATE VARIABLE 

THETI 

REAL 

INCIDENCE ANGLE FOP SURFACE REFlECTlON 

THTlLO 

REAL 

ERROR IN MEASUREMENT OF PHYSICAI ANGLE 

title 

character 

■ HEADINGS FOR PRINtOUT 

TWOPl 

REAL 


wavlth 

REAL 

WAVELENTH 

XI 

PEAL 

START OF SURFACE 

X2 

REAL 

END OF SURFACE 

XAC 

REAL 

POSITION OF AIRCRAFT 

XAO 

REAL 

STARTING POINT OF FLIGHT PATH SFGMENT 

XACF 

REAL 

ENDING POINT OF FLIGHT PATH SEGMENT 

XAl 

REAL 

POSITION OF image antenna 

XC 

REAL 

INTERCEPT OF SURFACE WITH LINE FROM AIRCRAFT TO ANTENNA 

XI 

real 

intercept of surface with NORMAL FROM AnTENNA 

XPEF 

REAL 

REFERENCE POlNT(NORMALLYsO.O) 

XS 

REAL 

SPECULAR POINT ON SURFACE 

YAC 

REAL 

POSITION OF AIRCRAFT 

YACl 

REAL 

STARTING POSITION OF FLIGHT PATH SEGMENT 

YACF 

PEAL 

ENDING POSITION OF FLIGHT PATH SEGMENT 

7.1 

REAL 

START OF SURFACE 

Z2 

REAL 

END OF SURFACE 

lt> 

REAL 

POSITION OF ANTENNA(X=0.0) 

ZAC 

REAL 

POSITION OF aircraft 

ZACI 

REAL 

START OF FLIGHT PATH SEGMENT 

ZACF 

REAL 

END OF FLIGHT PATH SEGMENT 

ZAI 

REAL 

POSITION OF antenna IMAGE 

ZI 

REAL 

INTERCEPT OF SURFACE WITH NOPMAl FROM ANTENNA 

ZRFF 

REAL 

REFERENCE POI NT (NORMALLY=0.0 ) 
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^UMR 


?54<)T* ACTI'/ITY * 


01 » RFPr-^T C''OF = 21* FECOPO COUM = 00010«> 


GLO^SAPY OF TFP“5;PP0GPA>'. “IA''6F 


AD 

RFAL 

SLOPE OF LINF ^OOM AIRCRAFT TO ANTENNA 

AIMAG 

REAL 


AR 

•REAL 

SLOP‘D OF LIMP FROM AIRCRAFT TO ANTENNA IMAGE 

RP 

real 

INTERCEPT OF LINE FROM AIRCRAFT TO A'.TEn:IA IMAGE 

COND 

RFAL 

ARRAY PSUFOOVY-' *^0P K.'PUT CONDITIONS 

CONDI 

RFAL 

ARRAY PSUFDONXW FOP INPUT CONDITIONS 

DELPHI 

REAL 

total p>j*ee frodr due to multipath 

DFLX 

REAL 

FLIGHT PATH INAPEMENT VALUE 

DFLY 

RFAL 

FLIGHT PATH INCRFMFHT VALUE 

DFLZ 

RFAL 

FLIGHT PATH INCREMENT VALUE 

FFOCUF 

real 

INDICATOR OP FDCUSED I NTERFECOmF TER ( 1 =F0CUSED) 

FXCT 

REAL 


I 

INTEGER 


INC 

INTFGER 

FLIGHT path INCREMENT LOOP CONTROL 

IFAIR 

INTEGER 

ANTF-MA PAIR 

IPOL 

INTEGER 

D0LARI 7/.T10N INDICATOR (1 = VERT . , 2 = HOR I Z . ) 

IRFFL 

INTEGER 

BINARY packed IMDICaTop OF WHICH SURFACES ARE REFLECTING 

I5FG 

INTFGER 

FLIGHT PATH SEGMENT LOO? CRNTROI 

I SURF 

INTFGFR 

REELECTING SURPACE LOOP COMTROl 

ITTY 

INTEGPR 

INDICATOR FOP PLOT SIZE 

J 

integer 


L 

integer 


LOGl 

LOGICAL 

LOGICAL INDICA'^OR TO DETERMINE IE REFLECTION EXISTS 

L0G2 

LOGICAL 

LOGICAL INDICATOR TO DETERNINE IE REFLECTION EXISTS 

LOG 3 

LOGICAL 

LOGICAL indicator TO DETERMINE IF REFLECTION EXISTS 

MCY 

INTEG'^R 


HFSl 

CHARACTF9 

PRINTOUT HEADINGS 

MFS2 

charactfo 

« M 

MFS3 

character 

ft M 

HFSA 

character 

ft ft 

WFS5 

charactfr 

ft It 

MINPH 

REAL 

INPHASE PORTION OF MULTIPATH 

HCUAC 

REAL 

CUADRATURP " • " 

NINC 

INTFGFR 

riuMRpR OF incremfnts this flight path segment 

NOUT 

INTEGER 

FILE CODE FOR LINE PRINTER 

NPAIR 

INTEGER 

number OF INTEPEEROMETER ANTENNA PAIRS 

NSFG 

INTEGER 

NUMPpR OF FLIGHt PATH SEGMENTS 

NSURF 

INTEGER 

MUMPER OF OPPLPCTING SURFACES 

OR 

BOOLEAN 

ARRAY CONTAINING RESULTS FOP A FLIGHT PATH SEGMENT 

our 

REAL 

ARRAY CONTAINING RESULTS FOP A FLIGHT PATH SEGMENT 

OUTMAX 

real 

MAXIMUM VALUE FOR CORRESPONDING ELEMENT OF OUT 

OUTMIN 

REAL 

MINI N'UM » » n . m 

PAID 

REAL 

ANTENNA PATTERN VALUE, IN DIR'ECTIDN OF AIRCRAFT 

PAIR 

REAL 

ANTENNA PATTERN VALUE IN DIRECTION OF REFLECTION 

PHA5 

RFAL 

PHASE OF GIVEN MULTIPATH AT GIVEN ANTENNA 

PHI 

REAL 

PHASE OF FRESNEL COEFFICIENT 

PHiRAR 

REAL 

AVERAGED ESTIMATE OF SINE OF PHYSICAL ANGLE 

phihat 

REAL 

ESTIMATE OF SINE OF PHYSICAL ANGLE FOR A GIVEN ANTENNA PAIR 

PPMTI V 

COMPLEX 

PERMITIVITY OP SURFACE 

0 

REAL 

DUMMY FLOATING VARIABLE FOR INPUT OF INTEGERS 

R 

COMPLEX 

FRESNEL REFLECTION COEFFICIENT 

RD 

DOUPLE 

RANGE FROM AIRCRAFT TO ANTENNA IN X-7 PLAME 

RD3 

DOUBLE 

RANGE FROM AIRCRAFT TO ANTPNNA 

RFAL 

RFAL 
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- 2I70T, /irTI'/tTY .Y : n. rorr r 21 , :?rc-f?r CDUrJ = 0002n9 

OOlOC pnrr~p/i>' ‘'AI'C,': 'r-jpuir-' lur. FPr-OPS P-jr TO SPrrU.AF PFcirCTlONS 

0020c AMP Frc'J?i‘;c. ffi A spArr-'-LLY avt'-acff r*)re’frpo‘'FTF9 for cp rp 

0030C 10 PAIPS OP «MTrMM4r , in PfFLFrTI'JP 5UPF '0^5 • AfO 1C PLIOHT PATH 

0032C SFOPFPTF 

0033 rOHMON/TAOF' /MOUT 

003A CnMN'ON/XfT/!TTV,-..rY»S<D(A) .F<CT(2) 

00 AO COMKOM 2* ( 10.2) »P( 10) «Yl ( 10) .21 ( 10) »Y2 ( 1C) <22 ( 1 C) ,PRP-,T I V ( 10) , 

00A56 FOH ( I'O) .XAf ! ( 10) ,YAC 1 ( 10) .7A0 1 ( 10) .YACF ( 10) .YACF no • ZAC F ( 1 0 ) . f' I NC ( 1 0 ) 

00A6 riMF'ISION COMPnO.Y) .COfPl (lO.A) 

00 A A POHJ 1/ALFACr (COPO.XAf 1 ) . (roe^o 1 ./) ) 

00A9 r lA'AMSlOK o>Jl'-“Y(lC) 

0050 r IMFNStOV SIIPF5 ( 10) .'-.UP FT ( 1 o ) « X A I ( 1 0 . 1 0 . 2 ) , 2 A ] ( 1 C . 1 0 « 2 ) , 

00556 FFLPHI <^l .OUF (1 00 . 1 5 ) .,0'lT m [ a (15) , V.'TK AX ( 15)' 

0056 P1KFN510N !PFFi_no) 

0060 Cny.PLFX PpvTiv.p(2) 

0062 rpAPACTFP P.F51-»6C/* A';Trf;'IA PAIPA*/. 

00636 PF52»60/' AroroAFT FLTOS'T PATH 5£0wfMT5*/. 

006A6 ►'F5.3»19/ 'OFFlFCT P'C-i SlIPF ACF 5 • / , 

00656 P'^5A'60/' ''.iRFaCF 5TAPT1MC COn^Ol VA TFS FT;Pi-.jO f OOPpiNATFS’' / . 

00666 PF.AATAB/TOFPWITI VITY 30I.IGHMF SS • / 

0067 CHAPACTfp TITLF»20{7) /*X AT GTA9T*.*Y AT START*. *7 aT START*. 

00686 *Y AT FINISH*. *Y AT FINISH*. *7 IT f j a> j Sh • , • MUt'P.FP f)F I \.CR PmFNT S * / 

0070, data ViAVLTH.TW''PI/. O'. 5. 2831853/ 

0071 PPAC NINPH.WlUAn 

0075. rO'JRLE PPFCISION PO ( 10.2 ) .Sa;T>J.R.7ff .dr 

0076 COUPLE PPPCISICN RO3(10.2) 

0077 LOGICAL LOG 1 . LOC-2 .L0G3 

0078 CATA IPOL/W 

0079 CATA XRFF.2PCF/6. 0.0.0/ 

0089C 

0090C *«*»*1NPUT ANTHNNA PAPAHpTrPS 

OlOOC 

0105 PFAn(20.150)0.-FFOCUS 

0106 IPOL=0 

0110 CO 190 IPAIPrl . 10 

0120 WPITF(11.130) IPAIR • 

0130 130 format (* INPUT ANTENNA MIGHTS FOP A.MTFNNA PAIR NUHPFR *.I2) 

OlAO RfaC( 20. 150)ZA UPAIP.ll .7A (IPAIR.2) 

0150 150 FORMAT (5F1C. 5) 

0160 IF(?A<IPAIR.l) .NE.O..OR.ZA(IPAIR.2 ).ne.O.)GO TO 190 

0170 NPAIR=IPAIR-1 

0180 GO TO 230 

0190 190 CONTINUF 

0200C 

0210C »»*»»input flight path PAPAMETFPS 

0220C 

0230 230 no 300 ISE6=1.10 
02A0 WPITE(11.250) I5FG 

0250 250 format ('INPUT AIRCRAFT X AND 7 COORDINATES AT START AND *. 

02526 ‘FINISH OF SFC-MFNT *.I2.*.AND NUNOPP OF INCREMENTS THIS SFf,*) 

0260 READ(20.150) XACl ( I5EG) .ZACl (ISEO) .XACF( ISEG) .ZACF (ISEG) . 

02626 0 

026A N1NC(ISEG)=0 

0265 WRlTEdl. 266) ISEG 

0266 266 FORMAT(*INPUT ST. AND FINISH Y COORD I NATES « SEGMENT *. 1 2 ) 

0267 READ(20.150) YACl (ISEG) .YACF(ISEG) 

0270 IF(S1NC(ISEG) .NE.O) GO TO 300 

0280 NSEG=1SEG-1 

0290 GO TO 3A0 
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0300 300 CONTINUE 
0310C 

0320C »**»«INPUT RE^L^CTIMG SI'PFA'"'^ 

0330C 

0340 ,3«0 ro 4<>0 I3URF=1.1C 
0350 WPlTF(ll»3f'0) ISURF 

0360 360 FOCwaTI • input X AND Z CPOPflNATFS \T ENOS CF SURFACE ‘.IZ) 

0370 RFAP(20«150) X 1 ( l SURF) .2 1 ( I SUPF) . X2 ( I SU9F ) .ZZ 1 1 «^UPF) 

0380 1F(X1 (I5UPF) .NF.X2(ISU0F).0P.?1 (ISURF) ,NF. 72 ( I SUOF ) ) 

03026 CO TO 410 

0390 NSURF=ISURF-1 

0400 60 TO 4R0 

0410 410 WPITFUI *420) ISUPP 

0420 420 FORMAT! • INPUT CnypLFX PfPNITIVITy AND ROUCH^'^SS PF SURFACE •* 
04226 12) 

0430 read (20* I 50) ORNT I V ( ISUPF) *F6H( I Sup'll 

0440 440 CONTINUE 

0450C 

0460C •••••COMPUTE SUFFACF PAPA''FTFRS 

0470C 

0471 480 ro 610 ISURF=1 .nSURF 

0473 IF(X1 ( I5UPFI .NE.X2(ISURF) ) 60 TO 4«6 

0475 SURFS( I5URFI =1.F37 

0477 SUOFI ( ISUPF) =1 .r37 

0478 DO 484 I PA IR= 1 *'!PA I R 

0479 DO 483 L=l*2 

0480 XA I ( I SL'PF* I°A IP *L ) =?.4X1 ( ISUPP) 

0482 ?AI (ISuPF*IPAIR*L)iZa(ipaiP*L) 

0483 483 CONTINUE 

0484 484 CONTINUE 

0485 60 TO 610 

0486 486 IF(71 ( I5URF) .nf. 72 (ISURF) )C-P TO 490 

0487 SURFS < ISUPF) =0.0 

0488 60 TO 500 

0490 490 SUpFS ( ISURF) : (7i ( i SURF ) -72 ( 1 SURF ) ) / (XI ( I SURF) -X 2 ( I SJRF) ) 

0499 500 CONTINUE 

0500 SURF I ( ISURFI = 7 I ( I SURF) -SURFS ( ISURF) 4X1 ( I SURF) 

0510 TEyoiri,.suCFS(ISURF)«*2 

0520 DO 600 I?AIF = l**;P«ir 

0530 DO 590 L = 1 *2 

0540 TFMP2= (?A ( ipair,l)-5URFI (ISURF) ) /TFmpi 

0550 XI=5UPFS(IS'IRF)»TEMP? 

0560 7I=7A(IPAIR,L)-Tf.m?2 

0570 XAI (ISUPF*IPAIR*L)=2.»XI 

0580 7A I ( I SUPF* IPAIR *L ) =2.«Z I-7A ( IPAI0*L) 

0590 590 CONTINUF 

0600 600 continue 

0610 610 CONTINUF 

0612C ••*»*PIJTDUT RUN CONDITIONS 

0614 '-■RIT.EnC*6l6)''FSl*7A*MFS?,(TlTLE(n .(CCND(J*I) *J= 1*10) *1 = 1*7) * 

06156 NES?*NES4*yF5‘' * (J* (COUOl (J*I ) .1 = 1*4) *PRFTIV(J) *F 6H ( J ) . .) = 1 . 1 D ) 

0616 616 FOPNAT(A60/’ UPPFR ' * 1 0 (Fft,2 *4X) / » LDWFP • * 1 0 (F6 . 2 . 4 X) // 

06176 A6O/6'(A2O*10(E9.3*2X)/) . A2C . 1 0-( I 5 * 6X ) / /A20/A56.14X.A35/ 

06186 • NUNRFR* *6X*2 ( *X* * 14X. V7* . 14X) /n( 3X* 

06196 I2*7(2X*fi7.5)/, , 

062CC 

0630C •••••increment a/C POSITIOr: AND EOmPUTF RANGF To FpRERFNCF POINT 

0640C 

0650 DC 1530 ISF6=l*NSEG 

0660 OELX= (XAC 1 ( IS^G.) -XACF ( ISE6) ) / Cine ( 1 5FS) - I > 

0565 DElVz ( YAC l ( I SFQ ) _Y ACF ( I S^F- ) ) / (NINE ( ISFG)-! ) 

0670 0ELT = (XACl ( ISFQ) -7 ACF ( I5EE-) ) / ( ') I NO ( I 5FC- ) - 1) 

0672 Si'rAN = 0.0 

0574 SNS0=0.0 

0575 DO 678 IPAlRz] *io 

0676 . IREFL ( IPAIP) =0 
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067B 678 
0660 
0686 
0690 
0695 
0700 
0710 
0712 
0720C 
,a730C 
0740C 
0750 
0760 
0770 
0772 
0775 
0777 
0780C 
0790C 
0800C 
0810 
0820 
08306 
0860 
0850 
0860 
0870 
0880 
0881 
0883 
0885 
0890C 
0900C 
0910C 
0911 

0920 920 
0930 
0960 
0950 
09606 
09706 
0972 
0980 
0992C 
0896C 
0996C 
1000 1000 
1010 
1020 
1030 
.1031 
1C60 
1050 
1052 
1060 
1O70C 
lOBOC 
1090C 
1100 
1110 

1120 1120 
1130 

1160 1160 
1150C 
1160C 
1170C 


continue . 

DO 1670 INC = ltNrNCMSE6l 
PHIB*R=0.0 

XACsXMTl nSEfil-(INC-l»«DELX 
YAC*YAC1 (ISEG»-nNC-l)*0EU7 
ZACaZACl ( 1 SEG) -( INC-1 1 *0CL2 
RRFF=DSORT( (XAC-XREF)«*2*lZAC-ZaEF)*«2» 

PATOsPAT (0. 0.1. 570796-ZAC/RREF* IPOD 

•••••select ANTENNA PAIR AN3 COHPUTE DIRECT hAN6E . 

00 1200 IPAIPsl.NPAlR 
DO 1160 L=l. 2 

RDUPA1P.D=DS0RT1XAC»»2* IZAC-ZAIIPAIR.L) )**2> 
RD3nPArR.D=DS0RT(RD<IPAIR.D**2*YAC«*2) 

MOUAO=0.0 

NINPHsO.O 

•••••SELFCT SURFACE 
00 1120 1SURF=1.NSURF 

RR=OSORT( (XAC-XAI (ISUPF.IPAIR.D )*»2+(ZACr 
ZAI (ISURF.IPAIR.Ln»»2) 

AP=(ZAI (tSURF.|PAlR.L)-ZAC)/(XAI ( I SURF • IPAIR .L > -XAC) 

BP = ZAI ( 1 SURF. IPAIR. L) -AR*XA 1 (ISURF. IPAIR. L) 

XS=(SUP.FI USURF)-BR) /lAR-SURFStlSURF) ) 

A0= (ZAC-ZA { IPAIR.tn /XAC 

XCs ISURFI 1 1 SURFl-ZA < IRA IR .L ) > / (AD-SURFS USORFI ( 
iF<SURFSnSURP» ,LT.1.E37)G0 to 920 
XC=X1(I$U»FI 
XSsXKlSURF) 

•••••oftfrmine if reflection occurs 

CONtiNUE. 

LCG1 = X1 (I SURF). GT.X2< ISURF) 

L0G2=XS.GT.X2( ISURF) 

LOGSsXS.GT. XI USURP) 

IF( (LOG1.AND.LOG2.ANO..NOT.LOG3.0R..NOT.LOG1.ANO.-.NOT.LOGZ 
.AND.L0G3) .AND. ( (XC.GT.O. .OR. XC.LT.XAC) .AND.XAC.LT.0..OR. 

(XC.LE.0..0R.XC.GT.XAC).AND.XAC.GE.0.) ) GO TO 1000 
I9EFL(IPAIR)=0RUREFLUPAIP) .2** (9* I SURF) ) 

GO TO 1120 

•••••PINO REFLECTION COEFFICIENTS 

THETI=ATAN( (l.*AR«SURF5( ISURF) ) / ( AR-SURFS (I SURF) ) ) 
THETI=APS(ThETI) 

call CFRSNL (PRMTIVUSURF) .RGHUSURF) .THETI.RU) .R(2) ) 
RHO=CAPS(R(IPOL) ) 

IREFL(IPAIR)=0P(IREFL(IPAIR) .2*»(I5URF-1) ) 
PH1=ATAN2(AINAG(R(IP0L) ) .real (R t IPOD ) ) 
PHAS=PHI.(RR-RD(1PAIR.L) )*TW0P1/WAVLTH 
P»TP=PAT(SURF5 (ISURF) .theti.ipod 
RH0PsRH0*P0( IPAIR.L)»PATR/ (PATD*RR) 

•••••accumulate multipath errors 

HOUAO=KOUAO.RHOR»SIN (PHA5) 

MINPH=MINPH.RH0R«C0S(PHAS) 

continue 

DELPHI (L) =ATaN(M0UAD/1.*minPH) 

continue 

•ji^^^AVEPAGE phase' errors 
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1180 PHIH*TnP»lR) = {R03CIPAI<»»l)-RD3(IPAlR.2)-WAVLTH/T/JC!PT* 

11828^ lOFUPHI (D-OCLPHI (2)-) I /-<2A « I0AIR.2-)-ZA ( IPAIR,1 )-) 

118A IF(FF0CUS.Ea.l.)PHIHAT<IPAIR)=3nrHAT( tPAI9)*tI ,-,5»PHlHAT(rPAIR) 

11816 *(ZA(IPAlR,2)-2A(IPAIR.in/R!5(IPAIR,l) )*.5 *(ZA(IOa;ir,1)* 

11886 ZA(IPAIRt2))/P0(lPAIR.l) 

1190 PHIPAR=PHIPAR*PHIHAT<iPAIR) /NPAIR 

1200 1200 CONTINUE 

1220C 

1230C *»A««PH1PAR IS TME REOUCEO OHASE NFASUPENEnT SJCH THAT THE 

1260C ESTIMATE OF POSITION IS ARCSIM (PHIRAR ) 

1250C 

1260C •••••COMPUTE ANGULAR FRROR 

1280 SnTh=2AC/RREF 

1320 1320 ThTILO= (SNTH-PH l3AR)^SORTn.O-PHlMAR^*2) 

1330C 

13A0C »****ST0PE OATA fop FUTURE OUTPUT 

1350C 

1360 00 1380 J=5*NPAIP*A 

1370 0UT(INC.J)=PHIHATU-4) 

1380 1380 CONTINUE 

1390 0UT«TNC«1)=R9EF 

UOO OUT(INC«2)=ThTilO 

1610 0UTnNC«3) =SNTh 

1620 00T(INC.6)=PHIAAR 

1630 OUTmIN(1)=-T.E-37 

1660 0UTMIN(2)=-1.E-37 

1650 OUTMAX(l)=-l.E-37 

1660 OUTMAX(21=-l.E-37 

1662 SHEANs SMEAN*THTIL0/NINC(ISEG) 

1666 SMSO=SMSO*ThTIlD»*2/N1NC(I5EG) 

1670 1670 CONTINUE 
1680C 

1690C •••••OUTPUT results ThiS SES^ENT 

1»0C 

1505 SVARsSmS0-SMEAN«^2 

1510 WRIT£n0fl520) ISEG, ( (OUT ( l ,J) ,J:l ,15) . I = 1 f MNC / ISEfiJ ) 

1512 WR1TE(10«1516)SmEAN»SVAF.,IRFFL 

1516 1516 FORMAT!* MEAN EPRORs • .El 1 .6 «6X. 'VAO j anCE = • E 1 1 .6// 

15166 • multipath INOICATOPS'/ZZX. 10(07, 6X) ) 

1520 1520 .FORMATdHl//* RESULTS POP FLIGHT PATH SEGMENT 
15226 'RANGE SPACE ERROR SIN TP SIM E5T*.20X, 

15266 'UNAVERAGEO EST im»TeS» /SO ( IX .FR. 3. 2X ,E1 1 .6 , 1 3 ( I X .F7.6 )/ ) HI/) 

1530 WRITEnO.1532) 

1532 1532 FORMAT (IHl/* PLOT OF ANGULAR ERROR VERSUS R ANGE • /3X ,* ERROR (RADI ANS) • ) 
1536 NOUTslO 

1S36 iTTTsO 

1538 CALL 5TUF(0UT( 1,2) .OUT(l.l) .NINCdSEG) ) 

1560 CALL PRNPLT(0UTd,2) ,OUTd,l) ,0UTM)N(2) ,OUTwi»((n ,nJTMAX(2) , 

15626 OUTmAX(I) .NINCdSEG) ,1,1,50,2) 

1566 WRITE(10,1566) 

1566 1566 FORMAT(50X, 'RANGE') 

1568 1530CONTINUE 
1550 STOP 

1552 END 

1 560C 

1570C •••••SUBROUTINE TO COMPUTE FRESNEL COFFFICIFNT 

1580C 

1590 SUBROUTINF CFRSNL (PRM.RGH.Th.PV.RH) 

1600 complex RV.RH.RRm.AI 
1610 A1=C50RT(PP,M-5IN(TH)*^2) 

1620 RV= (PRmacoS(TH) -A1 ) / (PRMmCOS(TH) *61 ) 

1630 RH=(C05(TH)-A1) /(C0S(TH),A1) 

1660 RE0=EXP(-157.R137^(RGH*CPS(TH) )««2) 

1630 PV=RVARED 

1660 RH=RHARED 

1670 return 
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1680 END 
1690C 

1700C »«**»rtJNCTION TO COMPOTE ANTENNA OATTEPNS 

ITIOC 

1720 FUNCTION PAT(SL«TH<n 

1722 DIMENSION eW(2) /. 3.1. 2/«TILT(2) /.1»*0.0/ 

1730C 1() IS UPTILT of ANTrMNA 

1740 ALPH=1.570796-TH-ATAN(SL) *T1LT( I) 

1750C »**»«PATTEPN FOR OPTr^Uf* HORN 

1760 T=1.9»5IN(ALPH) /SINIRWd) ) 

0001770 PAT=SIN(T)/T 

0001780 IF(PAT.LT.0.)PAT=-PAT 

0001790 return 

0001800 END 

1850C 

1860C •***»MA<E more points TO PLOT 

1870C 

1880 subroutine STUF(Y.X.N) 

1890 DIMENSION X(l) .Y(l) 

1900 ISTUF=100/N 

1910 IFnSTUF.LE.llCO TO 1980 

1920 DO 1970 I=1.ISTUF-1 

1930 DO 1960 Jsl.N-1 

1940 Y(N*J* (N-D* (1-in = IY(J*1)-Y( J) ) /ISTUF«1*Y(J) 

19S0 X(N*J*(N-l)*<I-in = <X(J.l)-X(J) )/tSTUF»I-»X<J) 

1960 1960 CONTINUE 

1970 1970CONTINUE 

IR75 N=(N-1I«ISTUF*1 

1976 WR1TE(U»1978)N.iSTUF 

1978 1978FORMAT<* PfflNTSs • . I4.2X. • I STUFs • . I« ) 

1980 1980 return 
1990 END 
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APPENDIX C 


(English Translation of Reference 6) 


1.0 INTRODUCTION 

Two transmitting antennas are symmetrically oriented on a platform 
which revolves about a fixed point, 0, as shown in Figure C-1. A 
receiver located at a point in apace, S, (at a range much larger 
than the transmitting antennas separation) receives the superimposed 
signals from the two antennas. The superimposed signals contain 
information (which when compared with a separate coherent reference 
signal) provide information on the azimuth and elevation angle of 
the point S relative to the transmitting antennas. By an optical 
processing scheme, a cross-correlation is performed between a 
transparent film with a transmittance function imposed on it and an 
illumination function amplitude modulated by the product of the two 
signals received at S. The optical cross-correlation produces luminous 
spots on a viewing screen such that the position of the spot in a two- 
dimensional vector polar coordinate system of angle and magnitude of 
the vector corresponds to azimuth and elevation angle respectively. 


2.0 CALCULATIONS 

2.1 Receiver Signals 

The antennas are located at points M, and M2 and the receiver at point S 
(See Figure C-1) such that separations d are: 


dl = SM]^ 
d£ “ SM2 


(C-1) 

(C-2) 
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2.1 Receiver Signals (Cont'd) 

The signals transmitted at the antennas have radian frequency 
( = 2 F,F the frequency in Hertz) and phase difference, g . 

At S the two signals received are: 

A^exp{j(nt - fi 

A2exp{j(nt - n + e)} 

The product of these signals is: 

A A 

cos{2fit - ^ (dj+d^) + 0} + cos{^ (dj-d^) + 6l 


The low-pass filtered part of this product is: 

P(t) = A^cos{^(d^-d2) + 6} = A^cos{^dj-d2) + 3) 
where: J\ = transmitter wavelength 

c = propagation velocity (speed of light) 

Calculating dj^ -dt 2 


(C-3) 


In the three-dimension reference frame located at the transmitting 
antenna, the point S is described by the angles 9 and \p and the 
distance d = OS (Figure C-2) , Then: 


=“TsM^ - 5^2 = 


(C-4) 


When the point S is far from the transmitting antennas (d"^7R) 


2d Z dj^ + d2 


(C-5) 
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2.1 Receiver Signals (Cont'd) 


Thus ; 

d^-dl = (d^-d^ ) (dj^+d^) - 2d (C-6) 

And from equation C-4 we get; 

d|-d 2 = 2R cosi|<s cosQ (C-7) 

When the two antennas are turning around the point 0 with radian 
frequency w, the time dependent angle 0 can be reference such that: 

0 = wt + 0s (C-8) 

Equation C-3 then becomes 

p(t) = AqCOS{^ cos«j)S cos(wt+0s) + B) (C-9) 

2.2 Transparency Function 

Consider a disk with equally spaced, alternating bands of opaque 
and transparent material. (Figure C-3). Also, consider the disk^ 
turning around the center 0 with the same angular frequency w as 
the two antennas . 

To simplify the equations, assume that the transmissivity of the disk 
is sinusoidal (instead of square) such that along the x-axis the 
transmissivity function may be described as: 

0(x) = + cos (C-10) 

where ; 

o^is the angular position corresponding to the transparency 
band maximum opaqueness relative to the center of the disk. 
0(X) is always positive 
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Figure C-3. Transparency Function 
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2.2 Transparency Function (Cont'd) 

This disk turns with the same angular frequency w as the antennas 
and, therefore, does not introduce an erroneous phase difference 
which would result in erroneous ©g and ^ spot references. The 
orientation of the disk will be chosen with the bands perpendicular 
to the x-axis at t = o. 


Consider the point, M( f,G), when the disk turns past this point 
the transmittance function is: (for ^ _ pcos9(wt + 0^) 

L(t) = + cos cos(wt + 0 q) + a} 

3.0 STUDY OF THE CORRELATION (SIGNAL PRODUCT AND TRANSMITTANCE ) 
Preliminary: 

(Z) is the Bessel function of the zeroth order with the 
following relationships: 

P’^cob(z sin 0) d0 = i cos (z sin 0) d0 = ^ / cos(z sin 0) d0 = Jo(z) 

2 ''^ 0 


1 ;^ 

— 0 sin (z sin 0) d0 


i. ^ ^ ain (z sin 0) d0 


1 

sin(z sin 0) d0 

2TT 


= 0 


The disk is illuminated with a single source of light with amplitude, 
or intensity, P (t) proportional to the product of the two antenna 
signals. Biasing the original signal P(t) so that the illumination 
function is never negative, we obtain for Ap >0: 


Ao{l + cos C cos ips . cos (ujt + 0s) + BD ) 

A 


(c-ii) 


(C-12) 


(C-13) 
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3.0 STUDY OF THE CORRELATION (Signal Product and Transmittance) (Cont'd) 
On the other side of the disk, at point M( P, Oq) , we can observe 
the product of that function with the disk's transmissivity function. 
Integrating the two functions as the disk revolves (which is what the 
eye sees) generates the two dimensional cross-correlation function 


G(p 0o) = ^ 2 ®°. {1 + cosC-**^ cos tJ/8.C08(u)t+es) + 6] } • { l+co8[^os(a)t+0o)-Kx] } dt 
» T 0 A 0 


(C-14) 


On setting wt = 0- and T = -■ we get; 

w 

G(p*0o) = AoBc C 1 + Ii + l2 + I3 H 
a) Calculating 1-^ ; 

I = C080 C08 [ 3 ^ cos( 0 + 0 p )] d 0 -£ / sin a . 8in [ ^ co 8 ( 0 + 0 o) 3 d 0 

^ 2'f 

we get: 

II = COSO • Jo 


b) Calculating I2: 


I = 1 /\osB-cosC-^ cosi|;„-cos(8+0,) ] d0-_L / slne-sinC;^ cos ^s.cos( 0 + 0 g)] d8 

2 — — o AS 2 '” ^ 


2tr 


we get: 


I2 = cosB . Jo ( ^ cos ) 


c) Calculating I3: 

2 TT 


211 _ 

I3 = cos cosi|/s . cos(0+0g) + 6] • cos cos (0+0o) + a] d0 


2Tt 


is made up of four integrals: 


(C-15) 

(C-16) 

(C-17) 

(C-18) 

(C-19) 

(C-20) 


f • j 

I3 = cosa . cosB . I - COSO. sing. I2 - cosB.sina.I + sina.sinB.I 


(C-21) 


1 
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3.0 STUDY OF THE CORRELATION (Signal Product and Transmittance) (Cont'd) 


which are: 

1 1 = ^ cos cos (e+e^)]. cos cos i/i^-cos (e+e^)] de (C-22) 

r’ _ A .2n cos cos (9+0 )]. sin cos -cos (e+e )] de (C-23) 

*0 ^ — J 0 0 ASS 

^ 0 

lo = ^ sin cos (e+e )]. cos cos -cos (e+e.)] de CC-24) 

I4 = 77/^^ sin cos (e+e )]. sin cos i|/.-cos (e+e^)] de (C-25) 

^ 0 ^ ^ ASS 


For all the integrals we make a change of variables, 6- for 0 and: 
4 ttR .A i A' = ^ - ^ cos ij/^.cos (6g-ep) 


A = ^ + ^ cos l^^.cos (65-6^) 


B = ^ cos i(/g.sin (65-6^) 


B' 


6 X 

cos .sin (65-6^) 


(C-26) 


B * 


Calculating Ij ; 


2it-9 


M ^ cos [i^^+B^.cos (e-i|))] de + 2— ; co? [v'^'^+B'^ cos (e-t|)')] 


47 


de 


(C-27) 


we get: 


Calculating 

de - ^ [> 7 ^.cos (e-*')] 


(C-28) 


-e. 


we get: 


de 
(C-29) 


I,=0 


(C-30) 
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Calculating 


I M 2tt- 6 yz r /~? 5“ 

13 = -^/ sin [v/a +B . cos (e-i|()] de + ^ ; sTn [a/A''^+B'^. cos (e-i|<')] de 

“0 


2tt-6_ 
- 6 . 


we get; 


13=0 

Calculating 

2tT:-0, 


* 4 ‘ 3 lr^ CO? [■/^'^•B'^.cos(e-*')]de - 4 - / co? [✓A^tB^.cos(e-i 

-% ->o 

we get: 


ij<)]de 


^4=7 (/aW)] 


Substituting these back into equation (C-21) we obtain: 

^3 “ I” cos (a- 6 ). Jq(\/A'^+B'^) + ^ cos (a+ 8 ), (\/a^+B^). ' 

and finally, from equation C-15: 

6 (p, 6 o)=A^jBo[l+COSa.Jo(^)+COSB.Jo(^ COSiJ<^ )+COS ( a-B ) (v^ '^+B ' ^ )+cos (a+B ) . 

4.0 INTERPRETATION OF RESULTS ^ 

Equation C-36 shows that the cross-correlation function is composed of 
five parts: 

a . Interpretation of the first term 

The term, A^Bq, is a> constant (independent of all parameters) 
and indicates uniform disk illumination of the observer. 


(C-31) 

(C-32) 

(C-33) 

CC-34) 

:C-35) 

)] 

:-36) 
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b . Interpretation of the second term 



shows that the disk has a maximum illumination at 0. 


According to the values of p , there are other maximiams 
and, since this term is independent of the angle 0^, we obtain 
concentric circles around the center of the disk at 0. 

However, this term vanishes for a = tt /2. 


c . Interpretation of the third term 


A ■‘'L— cos 4> 


does not depend on P or 0 so it also 


shows a uniform illumination of the disk with the intensity 


depending on <|)g. If we set g equal to ir/2, this term goes to 


zero , 


V 


d. Interpretation of the fourth and fifth terms, 

They are : x 

^ 0^0 cos (a+6) Jg (\/a +B 


We set the vectors u and v (Figure C-4) with the values 

' ' /9 — 7 I 

A, A, B, and B to get \/PC+B = |u + vj 

- ->■ 

= |u - v| 

The vector v, which is fixed, describes the point corresponding 
to the receiver at S. u defines the observed point on the 
disk. 


(C-37) 


(C-38) 


The illumination resulting from the cross-correlation, G( P,0o)» 
then is a function of only u. This function is maximum for some 


P\Ja 


values, such as\jA^ + B^ = 0 
to Jo(0), 



which corresponds 
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However, we cannot dispose of both parameters a and 3 , so it is 
impossible to suppress all of the following: 

The concentric circles around 0 

The uniform illumination with intensity dependent on (f* s* 

The false spot 

We therefore have to make a choice of which to remove. 


CONCLUSIONS 

Retaining the false targets (since we cannot remove or suppress 
all three of the above factors) we will obtain an illuminated 
spot on the disk whose angular position is representative of the 
azimuth and whose position relative to the center of the circle 
(distance = 4 i5 R cos <f>) allows us to calculate (j)s, the elevation 

angle. 
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5.0 APODIZATION BY AMPLITUDE MODULATION 

We now study the case where the signals received at the point S 
are amplitude modulated by the rotation of the antennas located 
at and M 2 when we use horizontal dipoles (Figure C-5). 


The received signal modulation is of the form Aj^ cos (wt +4>) and 

equation C-3 becomes : p(t) = A [1 + cos 2 (wt + ij^)] cos cos ij/^.cos (wt + 6 ^) + B 

0 ASS 

The product of the transmittance function B(t) will have the same (C-39) 

terms as derived above in one part and the other part will have new 
terms originating from the development of: 

I=j,5^COS(2wt+2)i;){l-COs[^ COSt(/^.COS(wt+eg)+e]}.{l+COS[^ COs(wt+0^)+a]}dt 


We again find the same integrals as calculated above with the 
exception of the factor cos 2 (wt + + 0 s - 0 s) under the integral 

sign. All of these integrations which yielded the Bessel function 
Jq ^zj now yield a term cos 2 ( <1^+ 63 "^b) *J 2 jz.j The result 


(C-40) 


has the form: 


i [Jo(Zi)+cos2(»+ej-e^) j^|2.)] 


We can now play with the coefficient cos 2 ( i|^ + 0s "©o) vary 
the form of the spots on the disk; in particular, the first sidelobe 


which is the most damped for the value: 


COS 2 1 


Apodization Effects 


The width of the 


^o Z 


curve at the 3 dB point is 1.55 and the first 


sidelobe is 507o and the width of the 


J2[z] 


curve is 2.0 at 


the 3 dB point with a first sidelobe of 177o. Figure C -6 shows the 
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Figure C-4. Vector Relationship 



Figure C-5. Amplitude Modulation 

Through Antenna Rotation 
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behavior of the curves along with the classical Sin Z ambiguity 

Z 

^function. — The-ambigutty-functionnshtainea^is not cyclic in the 
optical plane; it is apodized in the direction 'P +Gs 

When it is transformed to the direction \p + 9^ -9^ = 12 it will equal 

Jp ^zj -J 2 j^zjwhich has increased secondary sidelobes. The phase of 
the amplitude modulation (which can be electronically generated, such 
as , at the receiver) with respect to the antenna rotation define s 
the direction of apodization in the optical plane. 
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